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EOlTlON  OF  1 NOV  A*  IS  OBSOLFTl 


This  report  is  the  first  issue  of  the  Air  Force  Academy  Aeronautics  Digest 


Future  plans  call  for  publication  of  the  Digest  twice  a year.  The  contributions 
represent  recent  work  reported  in  a concise,  scholarly  and  timely  manner  by  students 
and  faculty  of  the  Department  of  Aeronautics,  members  of  other  Academy  departments 
and  F.  J.  Seiler  Research  Laboratory,  researchers  directly  or  indirectly  connected 
with  USAFA-sponsored  projects  and  authors  in  fields  of  interest  to  the  USAFA.  Ir 
addition  to  complete  papers,  the  Digest  will  include,  when  appropriate,  abstracts  of 
lengthier  reports  and  articles  published  in  other  formats.  The  editors  will  consider 
for  publication  contributions  in  the  general  field  of  aeronautics,  including: 

-Aeronautical  Engineering 
--Flight  Mechanics 
--Propulsion 
--Structures 

-Fluid  Mechanics 

-Thermodynamics  and  Heat  Transfer 

Papers  on  other  topics  will  also  be  considered  on  an  individual  basis.  Contributions 
should  be  sent  to: 

Editor,  Aeronautics  Digest 

DFAN 

USAF  Academy,  CO  80840 

In  this  first  issue  it  seems  appropriate  to  mention  why  publication  of  this  digest 
was  started.  As  part  of  the  normal  Aeronautical  Engineering  curriculum,  the  Depart- 
ment of  Aeronautics  sponsors  a tremendous  number  of  student  projects.  Many  of  these 
projects  are  well  worth  permanent  documentation  for  possible  use  in  future  research 
programs.  We  also  feel  that  the  preparation  of  work  in  a concise,  journal  style  is 
conceptually  valuable  for  those  students  whose  work  is  worthy  of  publication.  No 
other  outlet  for  student  papers  in  aeronautics  now  exists. 

In  the  Aeronautics  Department,  and,  in  fact,  the  U.S.  Air  Force  Academy,  many 
faculty  members  are  involved  in  formal  research,  and  others  work  on  academic  problems 
for  use  in  the  classroom  or  for  personal  clarification  of  some  principle.  In  many 
cases  the  formal  researchers  have  acquired  minor,  secondary  or  preliminary  results  but 
have  not  come  to  the  point  of  preparing  a project  report.  In  the  case  of  those  invest- 
igators of  academic  problems,  many  of  their  studies  represent  original  work  not  prev- 
iously published.  While  the  extent  of  these  studies  may  not  warrant  submission  to  a 
professional  journal,  it  is,  nevertheless,  important  to  retain  a permanent  record  for 
future  use. 


I 

I n discussing  these  points  with  rapt  Donald  Itishop  ot  the  Department  ol  History, 
he  suggested  that  the  Aeronautics  Department  would  do  well  to  start  its  own  | ourna 1 
At  the  time  the  statement  seemed  humorous,  but,  as  time  passed,  the  idea  seemed  to  make 
more  and  more  sense.  The  result  ot  Don's  suggestion  you  now  see  he  lore  you.  We  hope 
that  this  Issue  will  he  the  first  ol  many  more  to  come. 

1’he  Aeronautics  Digest  is  edited  by  t apt  1 . 1.  lumper,  PhD  and  rapt  M.  I.  Power. 

PhD  with  the  cooperation  of  the  Department  of  Ingltsh.  in  p.n  t i on  l.i  i , the  Digest  would 
like  tvi  thank  tupl  k.  I.  Stevenson,  't.  who  served  as  liaison  ollivci  I I om  Dll  Ni  , gav 
many  helpful  suggestions  and  piovided  the  final  editorial  levtow. 
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Abstract 


Experimental  wind  tunnel  data  is  presented  for  a right  conical  blunt  nosed  model  at 
angles  of  attack  of  SO  , 135  , and  180  , through  a maeh  number  range  of  1.2  to  3.0.  Hie 
axial  and  normal  force  coefficients  obtained  from  the  tunnel  data  acquisition  system 
are  converted  to  lift  and  drag  coefficients  and  extends  the  theoretical  data  published 
by  the  National  Aeronautics  and  Space  Administration,  NASA,  beyond  the  zero  to  90" 
range.  The  results  of  this  paper  could  be  included  in  the  aerodynamic  routines  of  a 
ballistic  trajectory  reentry  vehicle  computer  code  for  more  accurate  trajectory  pre- 
dictions of  a tumbling  reentry  vehicle. 


1.  Introduction 

With  the  advent  of  Intercontinental  Ballistic  Missiles,  ICBM,  to  world  political 
confrontation,  defense  of  the  continental  United  States  becomes  tremendously  complex  and 
expands  the  frontiers  of  technology.  rtiis  paper  addresses  a small  portion  of  the  ICBM 
defense  problem. 

Consider  the  scenario  of  a first  strike  attack  by  Russia  on  the  United  States. 

Varly  warning  radar  will  detect  and  tract  the  Russian  ICBM1 s near  their  trajectory 
apogee, and  linked  with  high  speed  computers  with  trajectory  prediction  codes,  can  deter- 
mine point  of  impact  in  a matter  of  milliseconds.  The  known  ballistic  trajectory  can 
then  be  integrated  with  an  Antiballistic  Missile,  ABM,  guidance  control  system  to  launch 
and  Intercept  the  ICBM.  However,  it  is  likely  that  several  high  altitude  precursor 
nuclear  weapons  will  be  detonated  to  disrupt  cotrmuni  cat  ions  and  blind  radar.  Added  to 
this  ABM  detonations  will  further  clutter  the  atmospheric  war  zone. 

At  first  glance  these  air  bursts  should  not  bother  the  ABM  system,  since  the  tra- 
jectory of  the  ICBM  was  already  determined  by  our  early  warning  radars  and  trajectory 
predictions;  unfortunately,  associated  with  the  atmospheric  bursts  are  shock  waves  with 
pressure  and  density  discontinuities.  If  a time  staggered  ICBM  happens  to  intersect  the 
shock  waves  of  a previous  detonation  the  trajectory  will  he  significantly  altered. 

Hie  ICBM  reentry  vehicle  can  skip  off  the  high  pressure  behind  the  shock  wave  and 
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unintentionally  extend  Its  range  similar  to  the  NASA  fears  that  the  Apollo  Spacecraft 
returning  from  the  moon  might  skip  off  the  Karth's  dense  atmosphere  and  be  plummeted  into 
outer  space.  Alternatively,  at  the  correct  intersection  angle  and  with  the  increased 
drag  due  to  the  higher  density,  the  reentry  vehicle  will  fall  short  of  it3  desired 
target,  see  Figure  1. 


Figure  1.  Modi  tied  Ballistic  Trajectories  due  to  Kncounter  with  a 
Nucle.it  Shock  Wave 


This  signficantly  complicates  the  ABM's  guidance  control.  The  ICBM's  trajectory 
has  been  modified  ilike  a maneuverable  ICBMl  such  that  the  ABM  system  can  no  longer 
predict  an  intercept  point.  So  work  has  been  done  at  the  Air  Force  Weapons  l.aboratory, 
Kirtland  AFB,  NM  to  couple  a six-degree-of-freedom  trajectory  code  with  an  atmospheric 
nuclear  burst  mode l to  predict  modified  trajectories  for  reentry  vehicles  that  had 
Intersected  the  burst  shock  wave  (Ref.  11.  Good  results  were  obtained  In  predicting 
Circular  Error  Probabilities,  CEP,  for  reentry  vehicles  passing  through  the  Influence 
of  a nuclear  burst. 

A key  to  any  trajectory  code  is  the  vehicle  aerodynamics . Accurate  calculations  of 
lift,  drag,  and  pitching  moments  are  Imperative.  These  parameters  can  alter  the  re- 
entry angle  and  range.  The  lift,  drag,  and  moment  coefficients  can  be  determined 
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theoretically  or  experimentally  tor  a given  shape  as  a function  of  angle  of  attack,  Mach 
number  and  Reynolds  number. 

In  the  course  of  time,  position  stepping  a trajectory  code,  angles  of  attack,  Mach 
number  and  Reynolds  number  are  known  at  the  previous  time  interval  and  all  parameters 
can  be  then  extrapolated  to  the  next  time  interval.  Mach  and  Reynolds  effects  on  lift, 
drag,  and  moment  coefficients  can  be  scaled,  however  it  is  more  complicated  theoretically 
to  predict  how  these  aerodynamic  coefficients  vary  with  incidence.  It  is  well  documented 
how  lift,  drag,  and  moment  coefficients  vary  with  angle  of  attack  from  + 85  (Kef.  2). 

For  spinning  reentry  vehicles  the  gyroscopic  effect  limits  the  incidence  to  generally 
less  than  10  . For  the  case  where  the  reentry  vehicle  intersects  a nuclear  shock  wave, 

the  Incidence  will  always  be  less  than  mo  . Hr  us , for  the  trajectory  codes,  the  aero- 

dynamic characteristics  of  the  reentry  vehicle  can  be  calculated  and  the  ballistic  tra- 
jectory predicted. 

For  the  original  scen-trlo  proposed,  there  is  one  further  complication  that  this 
paper  specifically  addresses.  Many  of  the  Russian  1CBM  reentry  vehicles  are  not  gyro- 
scoplcully  spinning  and  have  a tendency  to  tumble  fevcn  a spinning  reentry  vehicle  could 
possibly  tumble  if  it  encounters  a blast  shock  wave  at  a critical  angle!.  Ill  us , in 
order  to  predict  the  aerodynamic  characteristics  of  Russian  reentry  vehicles  the  lift, 
drag,  and  moment  coefficients  must  be  determined  for  a full  3b0  angle  of  attack,  rather 
than  the  widely  published  > 85  angle  of  attack  data.  It  is  the  intent  of  this  paper  to 
relate  experimental  data  for  lift  and  drag  coefficients  for  a Russian  type  ICBM  reentry 
vehicle  at  angles  of  attack  greater  than  85 

11.  F.  xjio  rlmenta  1 Approach 

In  order  to  obtain  lift  and  drag  coefficients  from  mo  to  180  incidence,  solid 
aluminum  models  were  designed  and  these  were  constructed  by  the  Air  Force  Academy 
Instructional  Technology  Division.  Hie  model  dimensions  are  depicted  in  Figure  2.  The 
shapes  are  typical  of  the  Russians  SS-M  ICBM. 


1 * 

If 

fc 


Figure  2:  General  Model  Shape 
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Using  the  flat  plate  approximation  for  the  base,  the  literature  suggests  u drag  coeffi- 
cient of  1.60  (Ref.  3).  The  drag  was  taken  as  607.  of  the  maximum  allowable  axial  load 
of  the  balance,  which  Is  334  N.  Tlie  dynamic  reference  pressure,  \ p 'fi.  Is  established  by 
tunnel  operating  conditions  and  can  be  obtained  from  a chart  in  the  tunnel  operation 
manual  (Ref.  4)  for  the  worst  case,  which  was  91  KN/m“"  . Tills  determined  a model  base 
area  of  0.177  square  meters. 

Testing  was  conducted  at  four  Mach  numbers  for  each  of  the  three  angles  of  attack. 

The  Mach  numbers  tested  were  over  the  full  supersonic  range  of  the  Trisonic  Wind  Tunnel 
from  1.2  to  3.0.  The  typical  Mach  number  profile  of  a Russian  SS-9  reentry  vehicle  is 
Mach  20  at  100  KM  down  to  Mach  one  at  impact.  Thus  the  tests  covered  the  low  altitude 
reentry  phase.  Data  was  collected  with  the  assistance  of  Fred  Jayne  during  November 
1977.  The  balance  data  was  collected  and  automatically  reduced  by  the  Hewlett-Packard 
Data  Acquisition  System  In  conjunction  with  the  Burroughs  TC-500  computer. 

111.  Discussion  of  Results 


The  theoretical  NASA  data  Is  tabulated  for  a sharp-nosed  right  cone  at  Mach  9.6. 
The  NASA  data  thus  had  to  be  corrected  for  Mach  effects  and  for  cone  bluntness  to  be 
compared  with  our  model  data.  The  following  discussion  and  equations  correct  the  data 
for  Mach  and  geometric  effects. 


The  axial  and  normal  force  coefficients,  and  C^,  respectively,  for  the  cone  are 
plotted  at  Mach  9.6  and  a stagnation  pressure  coefficient  of  2.0  according  to  Newtonian 
Drag  Theory.  A better  approximation  for  the  stagnation  pressure  coefficient  can  be 
made  for  low  supersonic  Mach  numbers,  i.e^M  < 3 by  the  following  equation,  (Kef . 2): 


y + 3 
y + i 


(‘ 


where  y - 1.4  for  air 

M - free  stream  Mach  number 


2 

(y  + 3)  MZ 


) 


(2) 


To  modify  the  tabulated  NASA  axial  and  normal  force  coefficients  for  Mach  effects,  each 
number  is  multiplied  by  the  ratio: 


C 

p 

max 

2 

where  the  numerator  is  determined  from  equation  2. 


(3) 


A blunt  cone  can  be  generated  by  superimposing  the  axial  and  normal  force  coeff- 
lcents  of  a sharp  right  cone  and  a hemispherical  shape.  The  superposition  can  be 
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achieved  by  the  following  equations  (.  Kef . 21: 


CA  " CA  (l  - 

cone  \ 

0 

1 + CA  , 1 

aphe re 

HOC  H 

c 

(4) 

S ‘ 1 N ( 1 

cone  \ 

■ ) 

+ c:  ( 

N \ 

sphere 

sec*  H 

c 

(5) 

where  C,  - axial  force  coefficient  of  sharp  cone  and  hemisphere 
A 

0 - normal  force  coefficient  of  sharp  cone  and  hemisphere 

h,  a and  9^  are  Riven  by  figure  6. 


Figure  6:  Definition  of  li , a,  and  H 

once  the  axial  and  normal  coefficients  have  been  modified  to  reflect  the  blunt  cone 
geometry,  they  can  be  resolved  Into  lift  and  drag  coefficients,  l'j  and  r , respectively, 
by  the  following  equations: 

* C.  cos  Of  + 0„  sin  « lb! 

D A N 

C,  * - C,  sin  o d C„  cos  O 1 71 

1.  A N 


where  o - angle  of  attack. 

flic  wind  tunnel  results  ure  plotted  in  figures  7-13  along  with  the  corrected  theo- 
retical data  provided  by  NASA  (Kef.  2).  11  ic  theoretical  data  Is  plotted  for  angles  of 

attack  from  zero  to  83  and  then  the  experimental  data  extends  the  curves  from  l*(V  to 
ISO'  . Both  lift  and  drag  coefficient  data  are  presented. 

Hie  experimental  data  blends  in  quite  well  with  the  NASA  theoretical  data.  1'he 
experimental  data  shows  the  highest  drag  when  the  model  Is  pointed  rearward  Into  the 
flow  which  Is  what  one  would  expect.  Hie  drug  drops  somewhat  near  the  '10  Incidence 
which  is  also  plausible  due  to  a reduction  lit  base  drag. 
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I'hi"  lift  coefficient  reaches  a maximum  near  1 IS1  angle  of  attack.  This  Is  a result 
of  the  large  flat  plate  area  of  the  base  being  exposed  to  the  free  stream.  The  lift 
coefficient  is  expected  to  return  to  zero  at  180'  incidence  because  of  axial  symmetry J 
•he  experimental  data  confirms  this  requirement. 

IV.  Conclusions 

Although  only  three  incident  angles  at  four  Mach  numbers  were  tested,  general  con- 
clusions can  be  drawn  since  we  expect  the  experimental  portion  of  the  curves, plotted  in 
figures  7-11,  to  |oin  smoothly  with  NASA's  theoretical  data  and  to  be  drawn  smoothly 
between  points.  The  highest  values  of  lift  coefficients  occur  near  135  and  the  highest 
drag  coefficients  are  obtained  for  a rearward  orientation,  180''  incidence. 

The  data  contained  In  this  paper  can  be  used  to  predict  lift  and  drag  coefficients 
for  a tumbling  reentry  vehicle.  With  a proper  curve  fit  to  the  data,  the  results  can 
be  Incorporated  into  the  aerodynamic  subroutine  of  a ballistic  trajectory  reentry  code. 
More  accurate  results  can  then  be  obtained  by  including  the  reentry  vehicle's  orientation 
through  a full  IblV  . 

In  the  wav  ot  recommendations  for  further  studies,  more  data  points  should  be  tested 
•he  lull  range  of  angle  of  attack  from  zero  to  180,  along  with  many  more  Mach  numbers. 
Higher  Mach  number  testing  is  required  for  the  higher  altitude  portion  of  reentry  where 
the  Mach  number  Is  near  .'0. 
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Drag  reduction  of  typical  light  aircraft  is  discussed  in  this  paper.  The  principle 
means  of  drag  reduction  can  be  achieved  by  reducing  the  wing  area,  thus  reducing  para- 
site drag.  Reducing  wing  area  naturally  involves  other  modifications  if  the  aircraft  is 
expected  to  have  the  same  or  better  performance.  The  modification  constraints  and 
related  configuration  changes  are  addressed  here,  such  as  reduced  span  or  chord,  full 
span  Fowler  flaps,  and  the  use  of  the  recently  developed  general  aviation  airfoil  devel- 
oped by  Whitcomb. 

I . Introduction 

The  wing  typically  accounts  for  almost  half  of  the  wetted  area  of  today's  production 
light  airplanes  and  approximately  one-third  of  the  total  zero-lift  or  parasite  drag. 

Thus  the  wing  should  be  a primary  focal  point  of  any  attempts  to  reduce  drag  of  light 
aircraft  with  the  most  obvious  configuration  change  being  a reduction  in  wing  area. 

Other  possibilities  involve  changes  in  thickness,  planform,  and  airfoil  section. 

This  paper  will  briefly  discuss  the  effects  of  reducing  wing  area  of  typical  light 
airplanes,  constraints  involved,  and  related  configuration  changes  which  may  be  necessary. 

I I . Constraints  and  Benefits 

The  wing  area  of  current  light  airplanes  is  determined  primarily  by  stall  speed 
and/or  climb  performance  requirements.  Table  1 summarizes  the  resulting  wing  loading  for 
a representative  spectrum  of  single-engine  airplanes.  The  maximum  lift  coefficient  with 
full  flaps,  a constraint  on  wing  size,  is  also  listed.  Note  that  wing  loading  (at  maxi- 
mum gross  weight)  ranges  between  about  10  and  20  psf,  with  most  4-place  models  averaging 
between  13  and  17.  Maximum  lift  coefficient  with  full  flaps  ranges  from  1.49  to  2.15. 

Clearly  if  C.  can  be  Increased,  a corresponding  decrease  in  wing  area  can  be 
max 

permitted  with  no  change  in  stall  speed.  If  total  drag  is  not  increased  at  climb  speed, 
the  change  in  wing  area  will  not  adversely  affect  climb  performance  either  and  cruise 
drag  will  be  reduced. 


Though  not  related  to  drag,  it  is  worthy  of  conment  chat  the  range  of  wing  loading 
in  Table  1 tends  to  produce  a rather  uncomfortable  ride  in  turbulent  air,  as  every  light- 


WING  LOADING  AND  Civl,v 
^nAA 

FOR  TYPICAL  S INCLF- PNC  INI 

AIRCRAFT 

Aircraft 

W/S  PSF 

CL 

max 

Cessna  1 50 

10.2 

1 .73 

Cessna  172 

13.2 

2.15 

9 1 

U 

Cessna  182 

16.9 

2.03 

Cessna  210 

21 .7 

2.01 

v 

Beech  C23 

16.8 

1 .89 

4 

Beech  V35B 

18.8 

1.85 

i 

J 

Grumman  Tiger 

17.1 

1 .92 

*, 

Bellanca  300A 

20.6 

1 .64 

\ 

Mooney  M20E 

15.4 

1 .85 

•i 

Piper  PA-28- 140 

13.4 

1 .73 

.1 

Piper  PA-28- 180 

14.4 

1.51 

1 

Piper  PA-28-200R 

15.6 

1 .49 

i 

Piper  PA-32 

19.5 

1 .92 

WINC  LOADING  AND  CLmax 

FOR  TYPICAL  TWIN-ENGINF 

AIRCRAFT 

Aircraft 

W/S  - PSF 

CL 

max 

Beech  Boron 

25.6 

1.42 

Beech  Duke 

31 .8 

1 .64 

Beech  Queen  Air 

29.9 

1.78 

Cessna  310 

30.7 

2.02 

Cessna  402 

32.2 

2.02 

Cessna  421 

35.2 

1 .86 

Piper  Seneca  II 

21 .9 

1 .80 

Piper  Navaho  PA-31 -350 

30.6 

1 .66 

Piper  Navaho  PA-31  P-425 

34.1 

1 .93 
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piano  pilot  Is  wo  1 1 aware.  The  only  way  to  reduce  this  gust  sensitivity  Is  to  Increase 
wing  loading. 

Typically,  wing  loading  tends  to  Increase  as  performance  (cruise  speed!  Increases. 
This  Is  particularly  evident  In  Table  2 which  presents  data  for  twin-engine  aircraft. 
But  gust  response  Is  proportional  to  the  ratio  of  calibrated  cruise  speed  to  wing  load- 
ing , V^/tW/Sl]  and  thus  Improvements  In  ride  due  to  higher  wing  loading  are  partially, 

If  not  completely,  offset  by  higher  cruise  speed. 

It  Is  also  evident  In  Table  2 that  even  though  wing  loading  Is  higher  than  for 

single-engine  aircraft.  It  Is  translated  directly  Into  higher  stall  speeds.  Twin-engine 

high  lift  systems  produce  virtually  the  same  as  shown  In  Table  1 for  single-engine 

‘max 

airplanes.  Thus  there  appears  to  be  an  equal  potential  for  reduction  In  wing  area  of 

single  and  twin-engine  aircraft  by  employing  Improved  high  lift  systems.  How  to 

achieve  higher  Cj  for  light  aircraft  Is  discussed  later, 
max 

But,  assuming  for  a moment  th«t  improvements  In  are  available,  making  higher 

max 

wing  loading  possible  for  a given  alrpiOne  or  class  of  airplanes,  it  Is  important  to 
consider  how  the  wing  area  should  be  reduced.  The  easiest  and  most  tempting  way  is  by 
reducing  span.  Not  only  does  this  leave  the  inboard  wing  structure,  mechanisms,  and 
wing-body  Junction  unchanged,  but  It  reduces  wing  bending  moments  making  possible  a 
lighter  wing.  But  reducing  the  span  increases  the  span  loading,  thus  reductions  in 
parasite  drag  through  a decrease  in  wing  area  are  countered  by  an  Increase  In  induced 
drag. 


F ^ 


On  the  other  hand,  reducing  wing  area  by  a decrease  In  wing  chord  decreases  parasite 
drag  almost  In  direct  proportion  to  chord  decrease,  and  If  span  remains  constant  there 
Is  virtually  no  change  In  Induced  drag.  From  an  aerodynamic  point  of  view  this  is  most 
desirable,  but  It  Introduces  possible  structural  and  weight  problems  because  aspect  ratio 
increases  while  spar  thickness  and  internal  volume  decrease  if  the  same  airfoil  section 
Is  used. 


To  understand  the  potential  and  the  constraints  of  drag  reduction  through  wing  area 
reduction,  consider  the  following  simplified  analysis. 

Assuming  that  the  parasite  drag  coefficient  and  span  efficiency  factor  remain  un- 
changed, the  parasite  drag  Is  directly  proportional  to  wing  area  and  Induced  drag  Is 
Inversely  proportional  to  the  square  of  the  span.  then  the  wing  drag  at  any  given 
flight  condition  may  be  written  as: 


where  D„  is  the  reference  wing  profile  drag;  IX,  Is  the  total  drag  of  the  reference 
PK  WK 


wing.  The  span  and  chord  are  denoted  as  b and  c with  a subscript  N Indicating 
reference  values.  For  simplicity  an  untapered  wing  Is  assumed.  Normalizing  equation  (11 


with  respect  to  the  original  reference  wing  drag,  "y  • (fives 
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whe re  P 


V 


the  ratio  of  parasite  drag  to  total  drag. 


If  only  the  wing  chord  Is  reduced,  then  the  change  in  total  normalized  wing  drag  is 


dll  - P 


dc 
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Thus  the  percent  reduction  In  total  wing  drag  is  equal  to  the  percent  reduction  in 
chord  length  times  the  original  ratio  of  parasite  to  total  wing  drag.  Clearly,  the 
benefits  of  wing  area  reduction  increase  with  air  speed. 


Consider  a typical  light  airplane  with  the  following  characteristics: 


Cross  weight  » 2800  pounds 
Aspect  ratio  * 7.4 


Wing  area  - 174  ft* 


frag  coefficient  of  body  and  empennage,  l 
Wing  parasite  drag  coefficient,  C 


li 


tl 

0.004 


0.017 


°BVH 


3J 


Airplane  efficiency  factor,  e » 0.75 
Cruise  altitude  « 8,000  ft 


If  only  the  chord  is  reduced,  then,  as  shown  in  Reference  1,  the  resulting  normal- 


ized total  airplane  drag,  D , is  shown  in  Figure  1.  Although  substantial  drag  reduct- 


ions are  possible,  constraints  are  imposed  by  the  requirement  to  cruise  at  a reasonably 
low  ltft  coefficient  and  stall  margin,  and  to  keep  stall  speeds  low  enough  for  good  take- 
off and  landing  performance.  Even  with  these  constraints,  however,  significant  reduct- 
ion* in  wing  area,  cruise  drag,  and  gust  response  are  possible  for  today's  general 
aviation  fleet. 
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Co  analyze  the  of  luce  ol  reducing  span  while  holding  chord  constant,  differentiate 
equation  v 2 1 with  respect  to  span  b.  1'hen 
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P - 2 (1  - P>  — 
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db 
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dp 


Kor  decrease  in  span  to  result  In  a net  decrease  In  drag  the  condition  - — 

Ul> 


tor  b • b muse  be  satisfied. 
K 


This  Is  true  only  It 


In  other  words,  a reduction  In  drag  by  reducing  span  can  be  achieved  only  if  parasite 
disc  Is  more  than  double  the  Induced  drag  at  the  t light  condition  In  question.  While 
th.is  mav  be  satisfied  during  high  speed  cruise,  it  Is  rarely  true  during  a climb.  And 
when  P i a reduction  in  span  Increases  induced  drag  more  than  It  decreases  parasite 

drag.  lor  a capered  wing,  P must  be  even  larger  than  the  value  given  in  (.51  to  achieve 
drag  reduction. 


the  limit  to  favorable  span  reduction  Is  found  by  solving 
dll 

which  yields  0,  assuming  P 2 '.  Again  from  equation 

dB 

that 


for  the  value  of  — — 
\ 

(x'  It  Is  easily  shown 
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db 
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equation  (ol,  plotted  In  Figure  2,  establishes  the  boundary  ot  favorable  span  re- 
duction of  a constant  chord  wing  as  a function  of  the  reference  wing  parasite  drag 
ratio,  P. 


I',’.,  {technical  developments 

It  Is  clear  that  wing  area  reduction  can  be  achieved  only  It  corresponding  Increases 

In  (,'j  can  be  designed  Into  light  airplanes  in  a practical  manner.  Several  recent 
max 

developments  Indicate  rhat  this  Is  a very  real  possibility. 


One  promising  development  Is  a new  family  of  general  aviation  airfoil  sections. 
Two  members  of  the  family,  the  \ t W 1 - 1 and  dAtWl-2.  have  been  defined  at  this  1 1 mo . 


As  shown  in  Reference  2,  Che  characteristics  of  these  airfoils  are: 


high  Cj  compared  to  conventional  airfoils  (see  Figure  3) 

'max 

gentle  stall  characteristics 

fairly  thick  section.  The  CA(W)-1  is  177.  thick.  This  helps  to  maintain  spar 
depth  with  reduced  chord  lengths. 

very  little  increase  in  C^  at  climb  lift  coefficients  (see  Figure  4).  This 

o 

combined  with  decreased  wing  area  offers  the  potential  of  significant  Increase  in 
single-engine  climb  performance  of  twins. 

Another  interesting  development  is  the  recognition  of  the  efficiency  of  spoilers  for 

roll  control  on  light  airplanes.  Among  other  features,  spoilers  permit  the  use  of  full- 

span,  or  at  least  increased  span,  flaps.  This  will  increase  C^  with  no  change  in 

max 

airfoil  or  flap  geometry.  Several  light  airplanes  are  now  using  this  concept:  the 
advanced  technology  light  twin  (ATLIT),  a modified  Seneca;  the  Kedhawk,  a modified 
Cessna  Cardinal;  the  RSTOL  Seneca,  a modification  kit  developed  by  Robertson  Aircraft 
Corporation;  and  the  Mitsubishi  MU-2. 

Another  method  of  increasing  C is  to  increase  the  Fowler  action  of  conventional 

'max 

single-slotted  flaps.  This  can  be  done  with  very  little  increase  in  complexity  or 

weight.  Figure  5 shows  the  very  large  values  of  C^  (2-D)  which  can  be  obtained  with 

'max 

a 0A(W)-1  airfoil  using  a 307.  chord  single-slotted  Fowler  flap, 

IV.  Flight  Test  Results 

Additional  confirmation  of  the  ability  to  increase  Cj  through  both  airfoil  design 

'mux 

and  flap  design  has  been  demonstrated  in  the  Redhawk  and  ATLIT  programs. 

Table  3 , from  Reference  3,  shows  maximum  lift  coefficients  obtained  on  the  Red- 
hawk  by  using  a 307.  chord  single-slotted  Fowler  flap.  Note  that  the  flap  covers  only 
477.  of  the  wing  span. 

The  ATLIT,  using  full-span,  307,  chord  single-slotted  flaps,  and  a GA(W)-1  basic 

airfoil,  generated  the  high  lift  data  shown  in  Table  4 . Clearly,  significant 

Increases  in  C^  are  possible  for  this  class  of  airplane, 
max 

Finally,  Table  5,  shows  drag  data  generated  during  flight  test  of  the  Redhawk.  The 
most  significant  result  is  that  parasite  drag  was  reduced  10.57,  by  reducing  wing  area, 
thickness,  and  span.  Tills  is  a significant  reduction,  and  it  illustrates  in  flight 
that  a reduction  in  wing  area  can  be  an  effective  and  practical  means  of  reducing  drag. 
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Configuration 


Redhawk 


Cardinal 


V , mpls  CL 

mo  \ 

V , mpF 
s 

Guise 

79. A 

1 .40 

64.7 

Kruger  flaps  only 

69.8 

1 .82 

Fowler  flaps  10° 

71 .2 

1 .75 

Fowler  flaps  10° 
and  Kruger  flaps 

62.8 

2.25 

Fowler  flaps  40° 

(30°  for  Cardinal  1 

64.4 

2.14 

55.0 

Fowler  flops  40° 
and  Kruger  flaps 

56.0 

2.83 

Notes  1 . Gross  weight  2500  lls 

2.  Redhawk  c.g.  location  7.2%  m. a. c.  (109  in. 1 

3.  Cardinal  c .i) . location  !9%m.o.c.  (Ills'. 3 ;n>) 

.’  l 


AT1.1T  mi  lMINAM  SIAM  DATA 


r 

L f 

Vc  -MPH 

->0 

CL 

max 

0 

76 

1-81 

10° 

66 

2.40 

20° 

61  .5 

2 .77 

u> 

o 

o 

59.3 

2.98 

40° 

59.4 

2.97 

Gross  weight  4200  lb 
Aft  c .g.  location 


Table  S 

I-OMI’ARISON  or  i'K.V  ■ I TiAKAi'T  I K 1 S I I T'-  I'l  l l 1 KM  I XI  I'  I ROM  K1  lOltr  TFST 


°P 

It  w 

P 

Cardinal 

Cruise 

0.0267 

4.67 

Full  Flaps 

0 .0462 

8.08 

Redhawk 

Cruise 

0.0380 

4.18 

Full  Fowler  and 

Kruger  Flaps 

0.0788 

8.67 
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takeoff  ground  roll  performance  model 

Fred  It.  Porter  III' 


Abstract 


An  acceleration-durlng-Cakeof f-roll  model  is  developed  using  two  unknowns,  a takeoff 
drag  coefficient  and  a rolling  friction  coefficient.  The  unknown  coefficients  are  then 
estimated  using  linear  regression  on  acceleration  versus  velocity  data  obtained  from 
four  different  takeoffs  with  varying  thrust,  weight  and  environmental  conditions. 

Finally  the  model  is  used  to  predict  acceleration,  time  and  distance  as  a function  of 
ground  velocity  for  standard  thrust,  weight  and  environmental  conditions. 

I . Int roduc t ion 


Traditionally  Hie  standard  takeoff  performance  data  of  an  aircraft  has  been  evaluated 
by  applying  individual,  first-order  corrections  for  non-standard  power  setting,  weight, 
temperature,  pressure,  and  headwind.  tills  method  of  evaluation  provided  an  estimate  of 
performance  on  a sea  level  standard  day  for  a specified  power  setting  and  aircraft  con- 
figuration. If  performance  was  necessary  under  other  than  standard  conditions  the  same 
process  was  applied  and,  although  tedious,  produced  reasonable  results  (Ref.  11. 

By  choosing  a general  linear  model  for  takeoff  acceleration  as  a function  of  ground 
velocity,  we  can  simplify  the  traditional  procedure.  Unknown  drag  and  rolling  friction 
coefficients  can  be  determined  by  using  a least  squares  fit  on  test  data  for  several 
different  takeoff  conditions.  We  can  determine  the  unknowns  (drag  coefficient  and 
rolling  friction  coefficient)  in  a general  takeoff  model.  With  these  coefficients  the 
model  may  be  used,  within  the  limitations  of  the  assumpt tons,  to  determine  takeoff 
acceleration,  time,  and  distance  quickly  and  expliclty  for  any  set  of  conditions. 

1 1 . The  Linear  Model 


A.  Assumptions 

Tlic  mode  1 incorporates  several  assumptions  as  follows: 

1.  Tlie  aircraft  remains  at  a constant  lift  and  drag  coefficient  throughout  the 
ground  roll.  Constant  landing-gear  strut  extension  fixes  the  aircraft  attitude  and 
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2.  The  runway  slope  Is  zero.  (Test  data  used  In  this  paper  were  taken  on  a zero- 


however,  the  model  may  be  modified  to 


slope  runway  which  simplified  the  regression 


accommodate  runway  slope.) 


3.  Not  thrust  Is  parallel  to  the  velocity  and  remains  constant  throughout  the  take 


off  roll.  Any  dependence  of  thrust  on  velocity  was  assumed  to  be  slight  and  was  absorbed 


In  the  regression  coefficients 


Aircraft  weight  and  configuration  remain  constant  throughout  the  takeoff  roll 


5.  Wind,  ambient  pressure  and  temperature  remain  constant  throughout  takeoff  roll 


Any  side  force  component  is  negligible 


Under  the  assumptions  of  part  A,  the  significant  forces  affecting  an  aircraft  during 


takeoff  roll  are  shown  In  Figure  l, where  T Is  the  net  thrust,  L Is  the  lift,  D Is  the 


drag,  W Is  the  aircraft  weight,  u is  the  rolling  friction  coefficient,  and  x is  the 


d istance 


Summing  the  forces  horizontally  yields 


Using  the  usual  definitions  for  lift  and  drag,  Eqn  (l)  is  solved  for  acceleration. 


V,  yielding; 


v “ w ( T - cd  * P VSS  " m w + a cL  % P * s 


whore  g is  the  acceleration  of  gravity,  C is  the  drag  coefficient,  p is  the  air  density, 
V is  the  longitudinal  component  of  air  velocity  relative  to  the  aircraft,  S is  wing  area, 
and  is  the  lift  coefficient,  Eqn  (2)  is  rearranged  to  separate  the  unknown  coeffi- 
cients, and  the  parametric  coefficient  C ' is  substituted  for  - y Cj  • 

gp  V*S 

V“  W T - V -15“  ' « 8 (J) 

The  unknown  coefficients,  O'  and  y , can  be  determined  by  linear  regression  using  test 


The  longitudinal  component  of  true  air  velocity  is  the  sum  of  the  aircraft  ground 
speed,  V , and  headwind  component,  , of  the  wind  velocity.  Using  this  sum  Eqn  (3) 

becomes  : 


V = ST  - c 1 I V;’  t 2 V V t V3  I • „ e < 

W DS2W\G  GH  H / **  b v 

yielding  a statement  of  acceleration  as  a function  of  ground  velocity  and  headwind  for 
the  takeoff  roll. 


It I.  Regression  Coefficients 


The  resulting  coefficients  of  the  regression  process,  C ' and  y , are  determined  by 
fitting  one  or  more  sets  of  takeoff  data  composed  of  pairs  of  observed  accelerations  V, 


and  ground  velocities  V . Eqn  (4)  is  rearranged  as  follows  ; 

1 


T y , 0s 


S'  2*  ( V + 2 VG  VH  + V ) + M 


New  variables,  a^  and  h , are  assigned  for  simplification. 

• • 

I It 
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Once  the  unknown  coefficients  have  been  determined  by  linear  regression  as  described 
in  section  III,  Kqn  (4>  can  be  used  to  make  predictions  of  takeoff  performance.  A more 
convenient  form  of  Kqn  (41  for  prediction  purposes  Is  ; 


v*> s 


V * 


V *0  SVH 


VG  * 


s o 

— Vs  - + “ T 

2V  H W 


(10) 


10 


-I,'  V,'lfc.f 

Ilf  I 


Since  the  third,  fourth,  and  fifth  terms  remain  constant  for  a given  takeoff  prediction 
Eqn  (10)  can  be  written 


x‘  = aV  2 + bV  + c 
b b 


where  a,  b,  and  c represent  the  appropriate  constants. 


The  takeoff  time  t is  determined  by  integration  of  Eqn  (11).  Since 

..  ^ 

x dt 


H 


aV  2+bV  +c 
b b 


after  integration  (Ref.  3) 


b2  -4ac 


-i  r2cvb  i 

Lv^^acJ 


The  takeoff  distance  is  determined  in  a similar  manner.  Since 


= v 

dt  G 


fix'f 


aVg2  +bVG+c 


and  the  integrated  expression  is  : 


^ - xi  - [aV  + bVG  + C] 


‘ 2c  (ts  " 
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V . Comparison  of  Model  and  Data 


Velocity  and  acceleration  pairs  from  four  observed  takeoffs  are  shown  in  Figure  2. 
These  data  were  recorded  during  takeoff  tests  on  a Lockheed  C-140  Jetstar  during  a flight 
evaluation  conducted  at  Edwards  AFB  in  August  1961  (Ref,  4).  Different  combinations  of 
aircraft  gross  weight,  ambient  temperature,  ambient  pressure,  and  surface  wind  occurred 
for  each  test,  and  these  data  were  used  to  determine  the  best  regression  coefficients 
C 1 and  y . The  predictions  of  the  resulting  model  are  also  shown  in  Figure  2.  Details 
of  each  data  set  and  information  on  the  "goodness  of  fit"  are  contained  in  the  Appendix. 


Figure  2.  C 140  Takeoff  Acceleration 

J2 


Phe  uncertainty  of  the  average  takeoff  thrust  Is  clearly  evident  In  Figure  3.  The 
thrust-level  had  not  stabilized  on  data  set  2 until  nearly  lift-off  speed.  It  Is  quite 
possible  that  throttle  adjustments  were  made  after  break  release.  In  spite  of  the  un- 
certainties the  takeoff  roll  can  bo  predicted  quite  accurately*  Table  l cixnpares  the 
takeoff  roll  predicted  by  the  model  to  that  predicted  by  the  usual  data  reduction  method 
used  at  Edwards  AFB  (Kef.  O. 


Table  1 


Comparison  of  C l-tl  Takeoff  Performances 


Parameter 


takeoff  ground  roll* 
takeoff  time* 


Model  Prediction 


1340  feet 


Reference  -*  Prediction 


3870  feet 


31  secs 


Predictions  aie  based  on  sea  level  standard  day  conditions,  no  wind,  no  runway  slope 
an  aircraft  gross  weight  of  40,000  lbf,  a net  thrust  of  3,833  lb j pet  engine,  and  a lilt- 
of t true  ground  speed  of  115  knots. 


VI . Conclusions  and  Recommendat  Ions 


A reasonably  accurate  model  for  prediction  of  takeoff  acceleration,  time,  and 
distance  can  be  made  by  estimating  drag  coefficient  and  rolling  friction  coefficient 
from  observed  velocity  and  acceleration  data.  The  model  can  be  used  to  predict  perform- 
ance under  a variety  of  takeoff  conditions  with  equal  ease. 

Phis  method  of  modeling  should  be  evaluated  extensively  on  other  data  from  tests  on 
different  aircraft. 


The  following  applications  are  suggested: 


1.  Ttrls  model  should  be  Included  In  flight  handbooks  in  equation  form  for  direct 
computation  of  takeoff  performance. 

3.  Phe  model  should  be  programmed  Into  Inertial  measuring  systems  on  board  aircraft. 
The  takeoff  predictions  could  be  computed  In  a "take-off  mode"  during  takeoff  roll  to 
provide  the  pilot  continued  comparison  of  actual  performance  versus  estimated,  Decisions 


iJ 


co  abort  could  bo  made  early  significantly  reducing  the  accident  rat 


3.  The  model  should  be  Installed  as  In  recommendat ion  - to  reduce  fuel  consumption, 
thrust  used,  and  engine  wear  on  takeoff.  With  accurate  takeoff  predictions  displayed  to 
the  pilot  the  throttles  could  be  set  to  the  minimum  thrust  level  required  to  achieve  a 
safe  takeoff.  Tills  would  Increase  engine  life,  save  fuel,  and  reduce  the  risk  of  engine 
failure. 


Koto r one  or, 

1.  Lush,  Kenneth  J . , "S tandarlz at  ion  of  Take-off  Performance  Measurements  for  Airplanes," 
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WIND  TUNNEL  TESTING  OK  SCALED  HANG  GLIDERS 
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Abstract 


This  report  examines  the  feasibility  of  aqulrlng  accurate  aerodynamic  characteristics 
for  hang  gilders  using  scaled  models  in  a wind  tunnel.  First,  the  paper  discusses  the 
results  of  a literature  search,  next  the  parameters  and  equations  applicable  to  collect- 
ion and  interpretation  of  data  collected  In  the  wind  tunnel,  and  finally  the  results  of 
wind  tunnel  tests  conducted  In  the  II SAFA  J'  x J1  subsonic  wind  tunnel. 

I . Introduct Ion 


A.  Background 

There  Is  no  longer  a question  about  the  Important  role  the  wind  tunnel  plays  In 
aeronautic  design;  ever  since  the  Wright  brothers'  first  poweted  l light,  aerodynamicist  s 
have  appreciated  the  wind  tunnel's  contribution  to  understanding  the  aerodynamic  and 
control  characteristics  of  an  aircraft  design  prior  to  flight  testing.  It  comes  as  no 
surprise,  then,  that  there  are  those  who  have  wanted  to  apply  wind  tunnel  techniques  to 
the  design  of  hang  gliders.  What  does  come  as  a surprise,  however,  Is  the  presence  of 
an  Immense  inertia  .pposlng  the  use  of  wind  tunnels  in  hang  glider  design. 


At  least  one  manufacturer  of  hang  gliders  openly  expressed  a need  for  some  sort 
of  unmanned  testing  of  new  designs,  and  the  wind  tunnel  was  one  ol  the  options  explored. 
Hie  manufacturer  was  Sun  Sail  Corporation  of  Denver.  In  an  interview  with  Cl lder  Rider 
the  co-owners , Brian  Jensen  and  Scott  Rojohn,  comnonted  on  the  state  of  the  testing  art 
and  the  need  for  unmanned  testing  (Ref.  1).  Jensen  and  Rojohn  typified  commercial  hang 
glider  manufacturers;  they  were  not  well  versed  in  aerodynamics.  Untypical  was  the 
technical  expertise  of  Rojohn,  whose  background  was  in  electro-mechanics,  structures  and 
materials.  As  is  true  for  most  manufacturers,  both  men  had  years  of  hang  gliding  exper- 
ience and  were  well  versed  in  traditional  hang  glider  testing  - build  them  and  fly  them. 
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Captain,  USAF, 


Assistant  Professor  of  Aeronautics,  DFAN 


t'o  present,  handling  characteristics  came  from  so  called  "public  testing"  - II  tin' 
accident  record  was  low  that  meant  the  handling  characteristics  were  good.  In  the  words 
of  Jensen  this  meant  the  "burden  of  testing  ended  up  on  the  unsuspecting  public."  Jensen 
went  on  to  say  that  "old  designs  are  becoming  obsolete  before  they  are  really  understood. 
And,  because  of  that,  the  potential  problems  of  the  old  ones  are  being  transfered  to  the 
new  ones  and  added  to  all  the  new  things  that  are  coming  out  in  the  new  ones....  We  are 
going  to  have  a lot  of  bad  things  happen  If  we  don't  start  trying  to  understand  what's 
really  happening  in  the  design  of  these  things," 

Because  of  the  nonavailability  of  a wind  tunnel  for  testing  theli  designs,  Jensen 
and  Rojohn  designed  a mount  for  their  pick-up  truck.  Using  a fish  scale  they  obtained 
Information  about  the  forces  required  to  change  angle  of  attack  and  were  able  to  modify 
the  hang  gliders  bv  trial  and  error  to  eliminate  "bad  handling  characteristics".  litis 
rather  crude  approach  to  design  was  an  Important  first  step  and  Jensen's  Interest  did  not 
stop  there. 

We  first  contacted  Jensen  In  January  of  1977  while  trying  to  locate  a so  called 
"high  performance"  hang  glider  design  to  use  as  a model  for  wind  tunnel  testing.  Jensen 
was  extremely  cooperative  and  loaned  us  one  of  his  designs  called  the  Swl ft.  Jensen  has 
followed  our  work  to  see  if  our  results  warranted  enough  confidence  to  overcome  the  un- 
favorable inertia  concerning  tunnel  testing. 

li.  Scope 

Vhe  purpose  of  this  project  was  to  demonstrate  the  feasibility  of  using  the  wind 
tunnel  as  a tool  for  designing  hang  gilders.  The  approach  was  to  1)  locate  full  scale 
hang  gilder  data  from  any  available  source,  2 ) construct  and  wind  tunnel  test  scale 
models  of  full  scale  hang  gliders  for  which  data  was  available,  and  31  draw  conclusions 
based  on  comparisons  between  the  lull  scale  and  model  data. 

I l . Available  Data 

We  found  only  one  good  set  of  data  on  a full  scale  hang  glider-type  wing  configura- 
tion, NASA  Technical  Note  D-1946  (Kef . 2).  The  wing  type  was  a standard  conical  K.gallo 
wing  shape.  flu1  work  reported  an  entire  range  of  aerodynamic  data  taken  In  full  scale 


wind  tunnel  tests.  The  data  had  been  validated  through  actual  fllghl  tests  In  an  asso- 
ciate program,  thus  could  be  considered  actual  full  scale  flight  test  data. 
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Flight  test  data  for  the  swltt  was  not  nearly  as  complete  as  for  the  standard  R gallo 
wing.  We  were  forced  to  reduce  pertinent  data  points  from  rather  unsophisticated  and 
uncontrolled  observations  reported  on  In  advertisements  and  popularised  hang  glider 
art lc les . 


Also  refered  to  In  this  paper  is  work  reported  on  by  Robert  Lee  Bass  (Ref.  3)  who  ran 
wind  tunnel  tests  on  small  scaled  models  at  the  University  of  Texas  at  Austin. 


III.  Theory 


The  theory  used  In  this  project  was  rather  simple  basic  aerodynamics.  In  the  follow- 
ing we  will  simply  give  references  to  the  equations  used  for  reducing  the  swift  flight 
test  data.  Additionally  reference  will  be  made  to  part  of  the  data  reduction  technique 
for  reducing  our  wind  tunnel  data  simply  because  It  Is  not  readily  available  elsewhere. 


A.  Reduction  of  Swift  Flight  Test  Observations 

Reference  4 reported  the  maximum  glide  ratio  and  stall  speed  of  a slightly  more 
advanced  swift  than  we  tested.  In  order  to  get  a ball  park  estimate  of  the  actual  flight 


characteristics,  we  backed  out  the  approximate  maximum  lift  to  drag  ratio,  L , ^ , and 


the  maximum  coefficient  of  lift,  C using  the  following  equations  (Ref.  5), 

max 


71) 


max  glide  ratio 


<n 


and 


VoV.t* 


ill 


(.’1 


where  is  the  coefficient  of  drag,  W Is  the  weight  of  the  man/hang  glider  combination, 


O Is  the  air  density,  V ( ^ the  stall  speed  and  S Is  the  planform  area  of  the  actual 

swift. 
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B.  Wind  Tunnel  Test  Reduction  Parameters 


Of  concern  to  the  correct  reduction  of  the  wind  tunnel  data  and  the  Interpre- 
tation of  that  data,  which  will  follow  in  another  section,  are  the  calculation  of  the  mean 
aerodynamic  chord,  the  calculation  of  the  correct  pitching  moment,  and  the  concept  of  the 
Oswald  Efficiency. 

1.  Mean  Aerodynamic  Chord 

By  definition,  the  mean  aerodynamic  chord  (MAC)  Is  the  length  of  a line  drawn 
through  the  centroid  of  one  side  of  a wing  parallel  to  the  airplane  longitudinal  axis. 

In  the  case  of  a triangular  wing  (the  classic  Rogallo  shape  for  example)  the  centroid  of 
each  side  of  the  wing  is  simply  located  1/3  the  distance  from  the  root  to  the  wing  tip. 
Irregular  shapes,  however,  must  use  the  following  equation  for  determining  MAC 


where  c is  the  chord  at  location  y and  y is  the  distance  from  the  wing  root.  For  deter- 
mining the  MAC  on  the  swift  the  wing  was  divided  into  difference  sections  and  the  integral 


was  approximated  by  a sunmatlon. 


f C1 ^yl 

£ c.  Ay. 


where  c^  is  the  average  cord  for  section  of  width  £ y^  . 

C.  Pitching  Moment 

Pitching  moment  Information  from  the  wind  tunnel  tests  is  given  referenced  to  the 
electrical  center  of  the  force  balance  in  the  tunnel.  In  order  to  get  meaningful  infor- 
mation the  pitching  moment  curve  must  be  mathematically  shifted  to  the  location  on  the 
model  corresponding  to  the  center  of  gravity  of  the  prototype. 

Of  interesc  to  this  project  was  the  slope  of  the  pitching  moment  coefficient 
curve  when  pitching  moment  coefficient  was  plotted  against  the  coefficient  of  lift.  The 
coefficient  of  moment,  (T,  , is  defined  as 


« 


o 


C m — — — 

M VD  Vs  S c 

(.■hero  the  M wav  the  moment  at  a given  angle  of  attack.  C was  then  plotted  against  a 

the  same  angle  of  attack.  Having  C.  vs  C,  the  slope  In  the  llniear  portion  of  the  curve 

n L 

C could  be  found.  file  negative  of  C is  known  as  the  static  margin,  SM,  or 

'C  (' 

L Ll 


SM 


(6 


Generally  Is 


plotted  In  relation  to  the  center  of  gravity  of  the  prototype  in  these 


cases  stability  Judgements  can  be  made,  and  the  SM  has  a physical  significance  tit  Is  th 
distance  from  the  center  of  gravity  to  the  neutral  point,  the  point  where  the  moment  is 
zero,  divided  by  cl.  In  the  case  of  the  wind  tunnel  data,  however,  the  SM  obtained 
purely  from  the  data  Is  In  relation  to  the  electrical  center  of  the  force  balance,  and  1: 
therefore,  the  distance  of  the  neutral  point  from  the  electrical  center  of  the  balance 
divided  by  c . In  order  to  find  the  SM  about  the  CG  location  for  the  prototype,  the 
neutral  point  must  be  found  by 


NT 


SM 


EC 


+ XKC 


(7 


Then  the  new  SM  for  the  CG  Is  found  by 


SMCG 


XNP  " XCG 


or  substituting  iqn  (71  into  : .jn  (81 


SMCG 


SMEC  c + XEC  - XCG 


(8 


(9 


We  were  then  able  to  use  the  SM  found  f rom  Iqn  (l,l  to  make  stability  Judgements  about 


I 


I).  Oswald  Efficiency 


If  the  coefficient  of  drag  for  a given  wing  is  plotted  against  the  coefficient 
of  lift  squared  the  linear  portion  of  the  curve  for  an  ideally  load  distributed  wing  is 
given  by 


C 


D 


1 


(10) 


where  .-V  is  the  aspect  ratio  of  the  wing.  For  a non-ideally  loaded  wing  Kqn  (10)  is 
modified  by  the  inclusion  of  the  Oswald  efficiency  factor,  e,  which  is  a measure  of  the 
efficiency  of  the  wing  (i.e.,  how  close  it  comes  to  an  ideal  wing)  where  e = 1 would 
correspond  to  an  ideally  loaded  wing.  Thus  Eqn  (10)  is  modified  as  shown  in  Eqn  (11) 


e rr  .V 


(11) 


For  conventional  lowspeed  aircraft  the  efficiency  factor  has  a value  of  about  0.91  to 
0.77  (Ref  3). 


IV.  Experiment 


A.  Pata  and  Observations 

Pnta  was  taken  on  four  different  flexible  wing  models  in  the  USAFA  2'  x 3' 
subsonic  wind  tunnel.  The  models  were  sting  mounted  and  strain  gauge  sensors  sensed 
sting  oriented  axial,  normal  and  torsional  forces.  The  signals  from  these  sensors 
together  with  angle  of  attack  with  respect  to  the  tunnel  axis  and  air  flow  information 
were  reduced  to  yield  Cj  , and  (about  the  electrical  center)  vs  angle  of  attack. 

Che  experiments  were  repeated  at  various  airspeeds  ranging  from  33  ft/sec  to  62  ft/sec. 
In  addition  to  the  force  coefficent  data,  observations  were  made  about  the  condition  of 
the  wing  at  various  airspeeds  and  angles  of  attack. 

B.  Models 

Four  wind  tunnel  models  were  tested.  The  models  were  made  from  two  planform 
areas:  Swift  (or  high  performance  planform),  and  Flexi-Flier  (standard  Rogallo  plan- 
form).  Each  of  these  two  planforms  were  tested  in  two  geometric  configurations.  In  the 
case  of  both  the  Swift  and  the  Flexi-Flier  the  configuration  change  was  a change  in  the 
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lodging  edge  shapes  i sco  Figures  1 and  ' 1 . Acrylic  sprav  was  also  used  Co  modify  the 
porosity  of  the  wing  material  i in  the  figures  unspraved  is  referred  to  as  "porosity"  and 
sprayed  as  "no  porosity"). 

V.  Results 


Oie  results  of  testing  are  presented  in  various  conventional  formats  in 
figures  * through  7.  In  all  of  these  the  angle  of  attack  was  measured  from  the  keel, 
t’he  titles  explain  the  contents  ot  each  of  the  figures.  Figure  d presents  the  Flexi- 
Klier  data  in  a modified  manner  adjusting  the  angle  of  attack  to  be  measured  from  the 
top  of  the  billow,  , instead  of  the  keel. 

h.  Observations 

l.  Flexi-Flier 

a.  Flutter 

In  all  configurations  the  Flexi-Flier  fluttered  excessively  at  angles, 
it  . from  -10  to  + IS'.  Flutter  was  present  Co  a small  degree  at  all 
angles  of  attack.  Starting  from  + 15°  and  reducing  the  angle  of  attack, 
the  flutter  started  at  the  trailing  edge  of  the  wing,  as  the  angle  of 
attack  was  decreased  the  flutter  grew  in  amplitude  and  frequency  as  well  as 
enveloping  more  ot  the  wing.  The  wing  was  totally  enveloped  by  + 4’  and  as 
the  angle  continued  to  decrease  the  wing  billow  reversed  becoming  more  or 
loss  stable  at  -10'. 

While  the  trend  was  the  same  in  all  Flexi-Flier  configurations,  the 
worst  flutter  condition  was  in  configuration  A with  no  acrylic  spray.  When 
the  acrvllc  spray  was  applied,  the  material  became  stiffer.  this  reduced 
the  severity  of  the  flutter.  In  conf iguration  H nylon  threaded  cellophane 
packing  tape  was  added  to  the  trailing  edge,  further  reducing  the  severity 
of  the  flutter. 

b.  Airspeed  Variations 

Airspeed  variations  did  not  seem  to  affect  the  data.  All  data  taken  at 
different  airspeeds  on  the  same  configuration  overlayed  each  other  exactly. 
Visual  observation  ot  the  billow  shape  and  later  comparison  of  photographs 
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taken  during  each  run  showed  identical  wing  shapes  at  all  airspeeds.  Hie 
flutter  had  a lower  frequency  at  the  slower  airspeeds  than  at  the  higher 
airspeeds;  the  amplitude  and  severity  of  the  flutter,  however,  were  unchanged. 
2.  Swift 

Hie  Swift  configurations  Included  wing  battens.  Those  battens  held  a 
cambered  shape  in  the  wing  at  all  angles  of  attack  including  the  negative  angles. 
No  flutter  occurred  with  the  Swift.  Data  and  wing  shape  remained  unchanged  with 
changing  airspeed . 

VI.  Discussion 

A.  Flexi -Flier 

1.  Coefficient  of  1,1  ft 

As  mentioned  in  Section  11  the  only  available  full  scale  data  for  comparision 
was  found  in  Reference  2.  Because  of  the  added  structure  in  the  NASA  test 
(manned  control  carriage,  etc.)  a comparison  of  the  drag  data  would  not  he 
meaningful.  Lift  data,  however,  was  compared.  Figure  i shows  an  overlay  of 
the  Flexi-Flier  configuration  data  and  the  NASA  data.  Some  simularities  are 
obvious.  The  slopes,  for  example,  are  nearly  identical.  There  is  also  an 
obvious  difference;  the  (1  equals  zero  intersects  differ  by  up  to  S degrees. 

It  we  examine  the  configuration  differences  we  notice  something  quite  inter- 
esting. the  main  difference  between  the  configurations  is  their  billow  heights. 

If  we  further  examine  the  profiles  of  configuration  A,  B and  the  full  scale 
wing  we  see  that  the  billow  angles  are  7 degrees,  1 degrees  and  1 ' degrees 
respectively  (See  figure  l).  As  mentioned  by  Bass  (ref.  3)  the  Cj  zero  inter- 
sect shifts  to  the  right  as  the  sail  billow  is  Increased.  Because  of  this 
statement,  we  decided  to  try  replotting  the  data  based  on  the  angle  between  the 
top  of  the  sail  billow  and  the  relative  wind,  rather  than  angle  based  on 

the  keel  line,  <1  . The  result  of  this  data  manipulation  is  shown  on  figure  S. 

Amazingly  all  sets  of  data  collapse  onto  the  same  line. 

This  result  suggests  two  things.  First,  it  would  appear  that  for  the 
Rogallo  type  wing  the  angle-of-nttack  used  for  plotting  vs  ft  should  be 
measured  from  the  billow  line  and  not  the  keel  line  (as  has  been  traditionally 


'>1 
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Clio  case).  Second,  when  tlie  data  Is  plotted  In  this  manner,  it  seems  apparent 
that  the  scaled  wind  tunnel  data  compares  very  well  with  full  scale  data. 

2.  Coefficient  of  Moment 


An  observation  can  also  be  made  about  the  scale  model  C slope  as 

CL 

compared  to  the  full  scale  model.  Figure  4 shows  the  original  data  about  the 
electrical  center  ami  the  adjusted  slope  to  the  assumed  center  of  gravity 
over  laved  with  the  full  scale  wing.  The  adjusted  curve  was  made  based  on 
equation  4 of  section  111.  As  shown  in  Figure  4,  both  the  model  and  the  full 
scale  wing  exhibit  near  neutral  stability  and,  bused  on  this  curve,  exhibit 
almost  identical  stability  characteristics.  It  should  be  noted  that  the  choice 
of  CO  is  somewhat  arbitrary  and  was,  in  fact,  chosen  to  make  the  curves  in 
Figure  8 coincide.  The  CO  location  needed,  however,  appeared  to  be  very 
reasonable. 

H.  Swift 

l.  Comparison  with  available  data 

As  mentioned  in  section  11  very  little  data  is  available  on  the  performance 

of  a full  scale  Swift.  From  advertisements  of  stall  speed  and  glide  ratio 

(i.e.,  best  forward  distance  to  decent  distance!  it  was  possible  to  back  out 

an  approximate  (1,/D)  and  C,  (Kef.  4).  The  actual  Swift  we  tested  was  an 
max  L 

max 

earlier  model  than  that  advertised,  the  Swift  A;  however,  since  the  newer  Swift 

advertisements  were  the  only  information  aval  lab le,  the  Swift  A information  was 

used  for  comparison  purposes.  Hie  main  difference  between  the  model  we  tested 

and  the  Swift  A was  a longer  wing  spun  reflected  in  the  Swift  A aspect  ratio  of 

5.28  compared  to  our  model  of  aspect  ratio  4.76.  Hie  performance  of  the  Swift  A 

was  thus  expected  to  be  somewhat  better  than  that  predicted  by  our  wind  tunnel 

tests.  Also  from  heference  4,  the  stall  velocity  for  the  Swift  A was  stated  to  bo 

13  Mi’ll.  Assuming  a 150  pound  man  plus  the  40  pound  weight  of  the  wing,  an 

area  of  170  square  ft  and  sea  level  conditions  Kqn  (21  yielded  a 0.  of  1.28. 

max 


According  to  heference  4,  the  max  glide  ratio  was  7 to  1.  Using  I'qn  ill,  (L/D) 

max 


SWIFT  CONFIGURATION  A 


SWIFT  CONFIGURATION  li 


FLEX  I FLIER 


Coefficient  of  Drag  versus  CoeDt lo lent  of  Lift  SqunreX  for 
the  Swift  Configurations  and  for'*'  Flext-Ftler  Configuration 


Figure  9 
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A SIMPLE  MODEL  OF  TEMPERATURE  GRADIENT  METAMORPHISM 

•k 

James  R.  Williams 


Abstract 

Increasing  winter  mountain  travel  by  skiers,  snow  mobllcrs,  and  climbers  means  an  in- 
creasing exposure  to  avalanches.  Slab  avalanches,  the  most  important  source  of  winter 
hazard  in  the  mountains,  occur  when  the  stresses  In  the  snowpack  exceed  the  strengths  of 
the  anchors  holding  it  In  place.  The  major  anchor  Is  the  shear  strength  between  the 
snowpack  and  the  surface  beneath  it.  Any  layer  with  a weak  shear  strength  is  a potential 
sliding  surface  for  a slab  avalanche.  Temperature  gradient  metamorphism  is  one  of  the 
primary  mechanisms  by  which  weak  layers  are  produced  in  the  snowpack.  A simple  model  is 
proposed  which  explains  the  mechanism  of  temperature  gradient  metamorphism  and  gives 
qualitative  results  which  agree  with  observations . 


I.  Introduct  ion 


Avalanches  seldom  touch  man  or  his  works,  but  when  they  do,  the  results  can  be  dis- 
astrous. From  data  reported  to  the  U.S.  Forest  Service's  Rocky  Mountain  Forest  and  Range 
Experiment  Station  it  has  been  estimated  that  there  are  about  10,000  avalanches  in  the 
United  States  each  year.  Typically,  l percent  of  these  harm  man  or  his  property,  account- 
ing for  an  average  ol  seven  lives  and  $300,000  in  property  damage  each  year  (Ref.  1). 

This  paper  will  discuss  the  two  major  types  of  avalanches  and  present  a simple  analysis 
of  snow  structure  metamorphisms , the  principle  factor  in  determining  the  probability  of 
slab  avalanches. 


It.  Types  of  Ava 1 anches 


Two  principle  types  of  avalanches  are  recognized:  loose  snow  avalanches  and  slab 
avalanches.  The  classification  is  always  based  on  the  snow  at  the  point  of  origin  with 
the  distinction  being  the  mechanical  character  of  the  snow. 


A.  Loose  Snow  Avalanches 


Loose  snow  avalanches  start  at  a point  or  over  a small  area.  Since  they  have 
little  Internal  cohesion,  loose  snow  slides  tend  to  move  as  a formless  mass,  growing  in 


size  as  they  descend.  A clearly  distinguishable  sliding  layer  is  not  always  present  and 
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the  line  separating  tin  sliding  snow  and  that  which  has  remained  in  place  is  often  in- 
definite. 

Loose  snow  avalanches  occur  when  snow  accumulates  en  slopes  of  steeper  angle  than 
its  natural  angle  repo  e.  Hiis  can  occur  by:  (1)  deposition  of  light  fluffy  snow 
under  conditions  of  little  wind;  ^ -”i  reduction  of  internal  cohesion  among  the  snow 
crystals  by  metamorphic  change-;  or  (.j'l  lubrication  from  percolating  melt  water.  Any 
small  disturbance  is  propagated  downward  from  crystal  to  crystal,  causing  the  snow  to 
slump  downhill  towards  its  natural  angle  ot  repose.  If  sufficient  momentum  is  attained, 
the  snow  may  slide  onto  slopes  of  lesser  angle. 

B.  Slab  Avalanches 

The  slab  avalanche  is  the  major  concern  because  it  forms  the  most  important  source 
of  winter  hazard  in  the  mountain--  . Slab  avalanches  are  characterised  by  the  internal  co- 
hesion of  the  now  Involved.  The  common  characteristic  is  a large  area  of  snow  begins  to 
slide  at  once.  A well-defined  fracture  line  where  the  moving  snow  blanket  breaks  away 
from  the  stable  snow  is  the  positive  identifying  characteristic.  A definite  sliding 
layer  is  usually  distinguishable. 


anchors.  At  the  beginning  of  a slide,  the  following  failures  typically  occur  (Fig.  1). 
Tl>e  crown  line  fails  in  tension  as  the  snow  is  literally  pulled  apart.  The  flanks,  bed 
and  stauchwall  fail  in  shear.  Since  the  bed  is  by  far  the  largest  area  involved,  it 


plays  the  major  role  in  holding  the  snowpack  on  the  mountain.  Any  layer  in  the  snow 
which  has  low  shear  strength  must  be  viewed  as  a possible  bed  for  a slab  avalanche.  In 
some  cases,  the  entire  snowpack  will  slide,  using  the  ground  as  the  bed  surface. 

Since  mechanical  forces  are  involved  the  structure  of  the  snow  is  important.  The 
snow  structure,  however,  does  not  remain  constant;  but  rather,  undergoes  metamorphisis 
which  greatly  effect  the  probabilities  of  slab  avalanches. 


III.  Types  of  Matt  amor  phis  is 


There  are  three  ^^pes  of  metamorphic  changes  which  alter  the  structure  of  the  snow 
crystals  after  they  are\in  the  snow  pack:  melt-freeze  metamorphism,  equitemperati^e 
metamorphism  and  temperature  gradient  metamorphism.  \ 


A.  Melt-Freeze  MetamorphXsm 


Melt-freeze  metamorphism ^curs  where  temperature  fluctuations  cause  melting  and  \ 
refreezing,  building  clusters  of  la^Ke  grains.  The  strength  of  a layer  of  melt-freeze 
grains  varies  widely.  During  the  meltAuart  of  the  cycle,  the  grains  are  essentially 
separate  and  are  held  together  only  by  siVface  tension.  Tlie  structure  is  then  extremely 
weak  and  ductile.  On  the  other  hand,  when  Jbe  layer  refreez.es  it  can  have  enormous 
Aren^th.  \ 


\ Eqtiltemperature  Metamorphism 


N^ot.weeh  temperatures  of  -TO°C  and  0°C,  ice  c r\\.i  1 s and  water  vapor  coexist  in 
the  s nowp.iOk,  amount  of  water  vapor  is  called  the  v^or  density  and  depends  not  only- 

on  the  t em p e r\t'a r e But  on  the  geometry  of  the  snow  crystal^^Convex  regions  have  a higher 
vapor  density  tmtri  conoawe  regions  (Fig.  2).  This  results  in  Abe  water  vapor  diffusing 
from  the  convex  areeA  of\he  crystal  to  the  concave  areas.  Additionally , the  water  vapor 


is  transported  to  the\:ianca\r>^  region  between  crystals  which  are  in  c^Mtact.  Fqultemper- 
ature  metamorphism  thus  \ends  ty  produce  rounded  grains  which  are  con n^^ id  to  each  other. 


rhis  results  la  a very  strong  snowpack. 


Convex 
high 
vapor 
dens  i tv 


Figure  - . Kquitemperature  Metamorphism  (Ref.  1) 
temperature  Gradient  Matamorphism 

tlie  temperature  dependence  of  the  vapor  density  (Fig.  implies  more  vapor  i- 
present  in  a warm  snowpack  titan  a cold  snowpack.  Since  vapor  tends  to  diffuse  from  area 
of  high  concentration  to  those  of  low  concentration,  a net  mass  transport  occurs  in  the 
snow  pack.  this  typically  results  in  mass  being  deposited  in  the  colder  layers  at  the 
expense  of  the  warm  layers.  The  enlarged  grains  that  result  are  characterised  by  flat, 
sometimes  steplike  or  straited  faces  which  intersect  at  sharp,  angular  corners.  Since 
there  is  little  cohesion  between  the  grains  the  layer  is  very  weak. 

IV.  Analysts  of  1'emperature  * radlent  Metamorphism 

For  simplicity  in  modeling  the  temperature  gradient  metamorphism  process  we  will 

o o 

assume  a linear  temperature  profile  with  0 C at  the  ground  and  -aO  C at  the  surface.  fh 
extension  to  arbitrary  temperature  profiles  is  straight  forward  ami  wall  be  explained 


t r o 

i 

Figure  3.  Variation  of  Vapor  Density  with  Temperature  (Ref.  2) 

Consider  in  detail  an  increment  of  snowpack  with  thickness  dx  (Fig.  4).  The  element 

* 

■ 1 

has  vapor  flowing  in  from  the  high  temperature  region,  vapor  flowing  out  into  the  low 
temperature  region,  vapor  being  produced  from  the  snow  crystals  and  vapor  being  condensed 
onto  the  snow  crystals,  A mass  balance  on  the  vapor  gives  ! 

■ 

M(x)  + M = M (x+dx)  + M (1) 

p c 

[s 

The  net  rate  of  crystal  growth  in  the  region  is: 


M - A = M (xl  - M (x+dx") 
c p 


O * density  of  the  mixture 
ptl,  0 + °AIR 


D 


diffusion  coefficient 


m ■ mass  of  water  vapor  per  mass  of  mixture 
°H,i' 

r — 

°ll.O  a °AIK 


Thus  it  follows  if  T(x)  > T(x+dx),  then  M(x)  > M(x-tdx),  From  eqn  (d)  we  see  this 
implies  a net  crystal  growth  in  the  area.  Since  crystal  growth  depends  on  M(x)  - M(x+dx), 
the  greatest  crystal  growth  will  occur  in  the  regions  where  the  vapor  density  gradient 
changes  most  rapidly.  For  our  linear  temperature  gradient,  this  will  be  the  region 
closest  to  the  ground  (Fig  51.  In  the  very  cold  regions  of  the  snowpack,  the  vapor 
density  gradient  approaches  zero,  so  little  temperature  gradient  metamorphism  will  occur. 
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Figure  5.  Mass  Concentration  of  Vapor  in  the 

Snowpack,  Linear  Temperature  Profile 


At  the  lower  boundary  of  the  snowpack  (i.e.,  the  ground), 

M(0)  = 0 


(4) 


Thus,  the  snowpack  consumes  itself  from  the  ground  up,  either  settling  or  creating  a 
hollow  space.  The  temperature  gradient  grain  growth  is  greatest  in  the  region  immediately 


adjacent  to  the  ground  and  decreases  with  increasing  distance  1 1 om  the  ground.  In  the 


very  c o Id  legions  ot  the  xnowpack,  there  is  essentially  no  temperature  gi  ah  lout  grain 
growth  at  all.  I’ll  l s corresponds  exactly  to  obsorvat  Ion. 

A particularly  interesting  case  arises  when 

M(x)*M(xtdxl  (S' 

Implying  Hi'  net  crystal  growth  in  the  region.  I'll l s Is  very  similar  to  the  case  s' t equl- 
t isnpei  at  me  metamorphism,  which  has  vapoi  being  produced  ami  consumes!,  but  with  no  net 
vapor  flow,  ts  It  possible  then  to  have  a non-uniform  temperature  sllstrlbut  Ion  and  on  i>  l - 
t omperat  m e met  amorphlsm  ’ 

For  no  net  crystal  growth,  Kqn  t1''  requires; 


Ktx  :?x ) ’Kix 


dm  I 
d x I 


xtdx 


An  approx  Imat  Ion  of  the  desired  temperature  piollto  (Klg.  b'  can  be  obtained  by  neglect- 
ing the  small  temperature  dependence  ol  and  considering 


dm  \ 
dx  )x 


dill 

dx 


) 


S td\ 


- constant 


If  this  temperature  profile  Is  present  lit  the  snowpack,  no  layer  will  show  a not 
crystal  growth  slnoo  the*  amount  ol  vapoi  diffusing  In  t vow  the  warm  side  will  match  the 
amount  of  vapoi  ill  (l  using  out  toward  tho  cold  side.  Of  course,  the  lavov  neai  (ho  ground 
will  still  ho  consumed  hv  this  process. 

It  Is  lutoiostfug  to  speculate  on  the  results  of  tho  mot  amoi  phlc  pi»voss  occui  tng 
with  this  particular  tomporaturo  profile.  Vho  physical  situation  can  ho  Imagined  as 
something  of  an  e^ul temperature  p rocoss  i s luce  thoio  Is  no  not  mass  gain  In  anv  rogioiO 
npon  wliich  Is  super  Imposed  a not  vapoi  flux. 


f'S 


would  simply  bo  reported  us  oqul t ompora t u ro  motumorphlsm 


Further,  It  seems  unlikely 


that  this  exact  profile  would  exist  frequently  In  nature. 


Since  there  Is  a vapor  flow  In  the  snowpaek,  the  second  possible  result  Is  temperature 
gradient  metamorphism,  even  though  there  Is  no  net  addition  of  mass  to  any  region.  This 
result  Implies  the  difference  between  temperature  gradient  metamorphism  and  equltempera- 
turc  metamorphism  Is  caused  by  the  vapor  flow  In  the  former  case  providing  condensation 
on  a preferred  side  of  the  snow  grain,  with  sublimation  occurlng  on  the  opposite  side 
(Fig.  71.  This  preferential  construction  would  result  In  temperature  gradient  grains, 
even  though  there  Is  no  net  mass  addition. 


Condensat Ion 


Vapor  Flow  (Constant) 

Figure  7 Preferential  Crain  Construction 
Caused  by  Vapor  Flow 

It  Is  known  that  temperature  gradient  grains  grow  In  a direction  opposite  to  the 
vapor  flow.  It  is  not  known  how  much  this  growth  depends  on  a net  mass  addition  to  the 
region  as  compared  to  how  much  it  depends  on  the  vapor  flow.  A series  of  experiments 
can  easily  be  Imagined : 

t.  beginning  with  a uniform  snowpaek.  Impose  ttie  temperature  profile  of  Figure  b 
and  observe  the  resulting  process.  Varying  the  snowpaek  thickness  will  vary  the  temper- 
ature gradients  and  the  amount  of  vapor  flow,  but  vapor  flow  will  remain  constant,  with 
no  net  mass  deposition  In  any  one  region.  Titus  grain  growth  could  bo  studied  as  a 
function  of  vapor  flow. 

Beginning  with  a uniform  snowpaek,  Impose  a linear  temperature  profile  and  ob- 
serve the  resulting  process.  Varying  the  snowpaek  thickness  will  vary  the  temperature 
gradient  and  the  amount  of  mass  deposition  in  the  different  regions  of  the  snowpaek. 
fit  us  grain  growth  could  be  studied  as  a function  of  mass  deposition. 


These  experiments  thus  would  focus  on  temperature  gradient  grain  growth  as  a function 


of  vapor  flow  and  mass  deposition.  Unfortunately,  they  would  be  extremely  difficult  to 
perform. 

V.  Extension  to  Arbitrary  Temperature  1’rollles 

For  an  arbitrary  temperature  profile,  Table  3 allows  a plot  of  m vs  location  in  the 
snowpack  (Fig.  4).  From  Eqn  (3),  we  see  that  areas  of  steepest  slope  correspond  to  areas 
of  greatest  vapor  transfer.  Areas  of  greatest  change  in  slope  correspond  to  areas  of 
greatest  crystal  growth. 

In  practice,  temperature  gradients  of  0.2°C/cm  to  0.3lC/cm  are  considered  to  be 
strong  gradients  which  will  produce  significant  temperature  gradient  metamorphism  (Kef.  11. 
Using  the  model,  we  may  incorporate  not  only  the  effect  of  temperature  gradient  but  also 
the  effect  of  temperature  Itself  on  temperature  gradient  metamorphism.  For  example,  a 
gradient  of  0.3°C/cm  at  -40°C  will  not  produce  significant  temperature  gradient  meta- 
morphism, but  a lesser  gradient  (say  9.2°C/cm')  at  -5^  C would  cause  significant  crystal 
growth.  Obviously,  more  work  is  necessary  before  it  would  be  possible  to  do  quantitative 
predictions  with  this  model. 

VI.  Conclusions 

This  simple  model  correctly  demonstrates  the  dependence  of  temperature  gradient 
metamorphism  on  both  the  temperature  and  temperature  gradient.  Warm  areas  in  the  snow- 
pack  are  shown  to  be  consumed  with  the  snow  turning  to  vapor  which  then  diffuses  to  the 
colder  areas  and  condenses.  This  can  either  cause  hollow  layers  or  settlement,  both  of 
which  are  observed  in  nature.  Essentially  undisturbed  layers  remain  at  very  cold  temp- 
eratures . 

The  model  suggests  the  possibility  of  a particular  non-uniform  temperature  profile 
which  might  behave  in  a manner  similar  to  equitemperature  metamorphism.  Apparently,  this 
has  not  been  observed  in  nature.  However,  without  specifically  looking  for  this  process, 
it  would  probably  have  been  reported  as  an  equitemperature  process. 


b7 


Experimentation  is  needed  to  determine  numerical  values  to  use  in  the  equations  so 
quantitative  predictions  could  be  made  and  to  determine  the  relative  dependence  of  the 


temperature  gradient  process  on  vapor  flow  and  mass  deposition 


VAPOR  PRESSURE  & VAPOR  DENSITY  DEPENDENCE  ON  TEMPERATURE  (Kef.  2 ) 


PROPERTIES  OF  AIK 


Calculated  from 


T(°C) 

PHB0*N/n’*) 

O,,  0<K*/^> 

-40 

12.9 

.011970 

-30 

38.  1 

.033979 

-20 

103.5 

.088605 

-10 

260.2 

.214270 

0 

610.8 

.48473 1 

i 


Table  3 

PROPERTIES  OF  THE  MIXTURE 


T(°C) 

til 

— 

°mlx(K8/nr'') 

-40 

.007839 

1.5271 

-30 

.022860 

1.4864 

-20 

.059703 

1.4841 

-10 

. 137664 

1.5565 

I 

0 

_1 

.272730 

1.7773 
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PRESSURE  DATA  FROM  THE  OPEN  POR\  FLIGHT  TES'l 
WITHOUT  APT  INPUT  WINIxfy 
* \ 

S.  E.  Icardi 


Abstract 

Ground  and  flight  tests  were  conducted  with  the  Airborne  User  laboratory  ^ demon- 
strate operation  of  the  Airborne  Pointing  and  Tracking  (APT I System  in  the  open  port- 
configuration  with  the  input  window  removed.  Preliminary  calculations  were  performed  to 
insure  that  a rapid  decompression  would  not  occur  when  the  APT  eyelid  was  opened  agu'inst 
a pressure  gradient  and  that  the  air  flow  through  the  APT  would  not  cause  damage  to  \ 
equipment  or  harm  personnel.  Hie  preliminary  estimates  proved  conservative.  The  meas- 
ured dynamic  loading  inside  the  APT  was  low  and  the  cabin  pressure  transients  were 
relatively  slow  and  did  not  degrade  flight  or  test  crew  performance. 


I . Background 


Flight  testing  has  been  conducted  with  the  Airborne  l.aser  Laboratory  ( A l . I . > to  confirm 
the  feasibility  of  operating  the  Airborne  Pointing  and  Tracking  System  (APT 1 without  an 
output  window  (.Ref.  1).  It  was  the  conclusion  of  the  testing  that  the  open  port  con- 
figuration appears  to  allow  better  beam  focusing  and  more  beam  intensity  on  target  than 
the  window  configuration.  These  promising  results  have  led  to  further  investigation  of 
the  open  port  as  a candidate  configuration  for  the  APT.  The  open  port  tests  were  con- 
ducted witli  a pressure  seal  between  the  aircraft  and  the  AIT.  To  be  a viable  config- 
uration, it  was  clear  that  the  open  port  had  to  be  tested  without  an  input  window.  For 
follow-on  tests  the  input  window  was  removed  from  the  bottom  of  the  APT  and  operation  of 
the  modified  APT  was  demonstrated  throughout  a restricted  ALL  flight  envelope.  This 
report  discusses  the  very  limited  aerodynamic  data  measured  during  ground  and  flight 
testing  of  the  follow-on  configuration. 

I 1 • Objective 

The  objectives  of  the  testing  program  were  to  demonstrate  operation  of  the  APT  in 
the  open  port  configuration  with  the  input  window  removed  and  to  conduct  beam  quality 
experiments.  The  beam  quality  experiments  were  aimed  at  detecting  beam  degradation 
caused  by  the  Increased  flow  through  the  beam  path  in  the  AIT.  The  particular  objectives 
of  the  aerodynamic  test  program  were  first  to  insure,  prior  to  flight,  that  the  air  flow 
through  the  APT  would  not  cause  damage  to  equipment  or  injury  to  peisonnel,  and  second, 
to  measure  the  environmental,  i.e.  pressure,  change  in  the  cabin  when  the  AIT  eyelid 
was  opened  against  a pressure  gradient. 
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III.  Method  of  Test 


Testing  of  the  open  port  APT  with  the  input  window  removed  was  conducted  in  two 
phases,  ground  testing  and  flight  testing.  The  ground  testing  phase  was  included  to 
insure  that  the  experiments  could  be  performed  safely  under  flight  conditions.  During 
ground  testing  the  aircraft  was  pressurized,  using  the  aircraft  pressurization  system, 
to  sea  level' standard  day  pressure,  14.7  psia,  approximately  2.5  psi  higher  than  the 
local  ambient  pressure.  The  APT  eyelid  was  then  opened  and  the  APT  parameters  were 
measured  and  recorded  using  the  existing  APT  and  aircraft  instrumentation  systems.  A 
pitot  static  system  was  installed  at  the  base  of  the  APT  to  measure  the  dynamic  pressure 
and  air  velocity  in  the  opening  created  by  removal  of  the  input  window,  Figure  1.  It  is 
assumed  that  this  opening  is  the  minimum  restriction  that  the  air  flow  will  pass  as  it 
leaves  the  cabin.  Therefore,  it  sees  the  maximum  air  velocity  and  dynamic  pressure. 
Cabin  pressure  was  recorded  using  the  cabin  pressure  instrumentation  installed  in  the 
aircraft.  Time  histories  of  the  cabin  pressure,  air  velocity,  and  dynamic  pressure  were 
made  and  compared  to  predicted  values  to  insure  the  decompression  of  the  cabin  would  not 
harm  personnel  or  the  APt  System. 

Flight  tests  were  conducted  to  gather  data  on  APT  performance  and  cabin  pressure 
transients  as  the  azimuth  angle  of  the  APT  was  changed  and  to  evaluate  beam  quality  as 
effected  by  the  increased  flow  through  the  beam  path.  The  only  pressure  data  recorded 
in  the  air  were  the  cabin  pressure  transients  as  the  APT  eyelid  was  opened  and  the  APT 
was  scanned  in  azimuth.  The  flight  testing  was  conducted  as  described  in  Ref.  2.  The 
inflight  measurements  of  cabin  pressure  versus  time  and  APT  azimuth  are  presented  in  the 
results  section  of  this  report. 

IV.  Theory 

Initial  theoretical  calculations  were  made  to  confirm  the  feasibility  of  such  test- 
ing. Preliminary  analysis  was  run  to  calculate  the  cabin  pressure  transient  as  the  eye- 
lid was  opened  on  the  ground  at  Kirtland  AFB  with  the  cabin  pressurized  to  sea  level 
standard  day  pressure.  Additional  calculations  were  then  run  under  what  were  considered 
be  the  worst  case  flight  conditions.  Plots  of  cabin  pressure  versus  time  were 
using  the  theory  of  conservation  of  mass  as  outlined  below. 

Time  rate  of  decrease  _ Net  rate  of  outflow  of  mass 
of  mass  in  volume  (V)  " across  the  boundary 

- 5^  f p dV  = f p (v  • N)dA  (1) 
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time  (s) 
density  of  air 


surface  normal  vector 


where 


surface  area  (ft*2  ) 


volume  (ftJ) 
velocity  vector  (ft/  ) 


for  steady  flow 


or  one  dimensionally 


constant 


where 


velocity  magnitude  (ft/  ) 


Applying  the  perfect  gas  relation  (p  = pRT)  and  the  Mach  relation 
m then  becomes  : 


For  adiabatic  perfect  gas,  > = 1.4 


Now  m may  be  calculated  by  knowing  the  static  pressure  (p),  Mach  number  (M) , area  of 
interest  (A),  and  total  temperature  (T  ). 


As  a conservative  estimate  we  can  assume  no  pressure  losses  through  the  APT.  The 
static  pressure  at  the  outlet  then  corresponds  to  the  static  pressure  outside  the  opening 
in  the  APr.  By  the  isentropic  relationship,  the  Mach  number  is  related  to  the  outside 
pressure  (p>  and  the  cabin  pressure  (p  ) by: 
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We  now  have  estimates  for  M and  p and  only  need  to  estimat\  the  cabin  air  temperature 

(T  i and  the  surface  area  of  the  outlet  (A)  to  solve  for  A.  Vlhis  is  an  instantaneous 

° \ . 

value  of  m and  must  be  adjusted  as  the  cabin  pressure  is  reduced.  If  we  assume  that  m 

is  constant  over  some  small  unit  of  time  (At  ),  we  can  estimate  the  change  in  pressure 
in  the  cabin  by  again  applying  conservation  of  mass. 

Mass  flowing  out  of  * Net  decrease  in  mass 
cabin  in  time  A t in  cabin 


m A t c A Mass  * A 0 (Vol) 
o 


\ 
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From  perfect  gas  (p 


I \ 


Rv  substitution  then 


if  i is  constant 
o 


or  the  change  in  cabin  pressure 


lit  A t Ki 


( fc  ' 

lbf  \ 

\ lbm 

where  m = instantaneous  mass  flow  rate  (lbm  1 

s 


K gas  constant  for  air  ( 

r “ cabin  air  temp  ( Rl 
o 

At  » time  interval  (si 
Vol  cabin  volume  (ft  ) 


l>  - P - A P 
o o o 

new 


■[(V-) 


with  M4*  now»  we  solve  for  a new  n,  using  Equation  3.  We  can  set  up  an  Iterative  proc- 
ess to  solve  for  P at  each  Ji:  which  results  In  a plot  of  cabin  pressure  (P  ^ versus 
o o 


L# 
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COEFFICIENT  OF  PRESSURE 


time.  Figure  2 is  a plot  of  the  initial  theoretical  analysis  of  cabin  pressure  versus 
time.  Initial  conditions  are  similar  to  those  that  were  set  in  the  ground  test.  In  this 
calculation  cabin  pressure  (P^)  was  assumed  to  be  initially  at  sea  level  pressure  (14.7 
psi),  outlet  static  pressure  (P)  at  5000  ft  pressure  altitude  (12.2  psi),  cabin  temper- 
ature (T  ) constant  at  70° F and  outlet  area  tA)  11.45  in*".  Test  point  1C  from  Kef.  2 at 
o 

turret  azimuth  angle  110  was  selected  as  the  flight  worst  case  since  it  presented  the 
largest  pressure  differential  between  cabin  pressure  and  outlet  static  pressure.  The 
cabin  pressure  (P  ) was  assumed  to  be  5000  ft  pressure  altitude  (12.2  psi),  the  outlet 
static  pressure  (P)  was  estimated  at  8.0  psi  using  pressure  coefficients  obtained  during 
flight  certification  of  the  open  port  configuration,  see  Figure  3,  cabin  temperature 

o 

(T  ) constant  at  70  and  the  outlet  area  (A)  11.45  in  . Figure  4 show  the  results  of 
these  calculations.  The  primary  purpose  of  the  calculations  was  to  show  that  even  with 
instantaneous  opening  of  the  APT  eyelid  under  worst  possible  flight  conditions  there  was 
no  chance  of  an  "explosive"  rapid  decompression. 

An  additional  concern  of  safety  personnel  was  the  flow  velocity  around  the  opening 
in  the  bottom  of  the  API'.  The  velocity  field  around  the  opening  was  calculated  using 
conservation  of  mass  and  assuming  perfect  gas,  steady  and  isentropic  flow. 


Let  us  assume  that 


m ■ p Av  1 constant 


and  that  we  have  a sink  located  next  to  a solid  wall. 


Figure  5.  Control  Volume  for  Mass  Flow  Calculations 

i . . • am-  a half  spherical  control  volume  as  shown  in  Figure  5 above.  fhe  mass 
•>  w*  lit  the  sink  is  denoted  by  m.  This  mass  must  from  conservation  of  mass  be 
. -m»«  flux  through  the  half  spherical  surface  area  around  the  sink 


» H .v:j 


where  p 


« density  of  air  (Ibm/  j) 

r ■ radius  from  sink  to  spherical  surface  area  (ftl 
v = velocity  normal  to  surface  (ft/s> 

1 t is  possible  to  neglect  compressibility,  since  the  velocity  at  a very  short  radius 
from  the  sink  is  well  below  Mach  0.5. 

The  density  then  becomes  constant  and  equal  to  the  density  of  the  air  in  the  cabin. 
The  initial  velocity  profile  can  now  be  calculated  as  a function  of  the  radius  of  the 
spherical  area, 


v 


m_ 

P 2 it  r2 


if  the  .itial  mass  flow  m and  the  cabin  air  density  p are  known. 


Assuming  perfect  gas  relation 
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P 


T 

r r 


KT 

o 


o 
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where  T *■  cabin  temperature  ( K) 

P^  * cabin  pressure  (lbf/.^g  ) 

Figure  6 shows  the  theoretical  initial  velocity  distribution  as  the  radius  ot  the 
spherical  surface  area  is  increased.  It  can  be  seen  that  the  velocity  drops  off  very 
quickly  and  is  less  than  10  fps  at  10  inches  from  the  sink.  The  velocity  field  around 
the  opening  does  not  appear  to  be  a significant  safety  hazard. 


V . Ground  and  Flight  Test  Results 

Tlie  ground  tests  were  conducted  by  pressurizing  the  aircraft  to  sea  level  standard 
day  pressure  (14.7  pslal  and  then  opening  the  APT  eyelid.  The  cabin  pressure  was 
recorded  using  the  aircraft  system.  Static  and  total  pressures  were  recorded  at  the 
opening  in  the  base  of  the  APT.  Figure  7 is  a plot  of  the  change  in  cabin  pressure  with 
time.  It  is  interesting  to  compare  this  measurement  with  the  theoretical  estimate 
(Figure  21  for  the  ground  test.  Hie  theoretical  estimate,  though  based  on  several 
rather  critical  assumptions,  proved  to  be  a conservative  calculation. 
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I'he  velocity  and  dynamic  pressure  calculations  were  made  from  the  pitot-static  data. 
Figure  8 is  a plot  of  the  change  in  velocity  with  time.  It  can  be  seen  that  the  velocity 

never  exceeds  200  fps  which  at  a cabin  altitude  of  5000  ft  corresponds  to  Mach  0.18.  In 

the  theory  section  we  showed  that  the  velocity  dissipates  very  quickly  as  the  distance  is 
increased  away  from  the  opening.  Figure  9 shows  the  change  in  dynamic  pressure  with  time. 
Hie  maximum  dynamic  pressure  was  measured  to  be  0.27  psi,  the  dynamic  pressure  was  well 
below  any  structural  limit  on  the  APT. 

The  pitot  static  system  was  removed  from  the  base  of  the  APT  for  the  flight  test, 
since  propagation  data  was  required  on  this  flight.  The  cabin  pressure  was  recorded 
during  opening  of  Che  APT  eyelid  and  also  through  an  azimuth  scan  of  the  turret.  Figure 
10  shows  the  cabin  pressure  response  versus  time  for  a typical  opening  of  the  APT  eye- 
lid. A comparison  of  the  theoretical  results  (Figure  4)  shows  a favorable  agreement  in 
the  data.  This  pressure  response  is  similar  to  what  might  be  seen  if  the  turret  was 
rotated  at  max  speed  from  an  azimuth  angle  of  zero  degrees  to  90  degrees.  Figure  11  is 
a plot  of  die  change  in  cabin  pressure  as  azimuth  angle  is  varied.  For  this  data, 

azimuth  angle  was  rotated  slowly  so  that  cabin  pressure  was  essentially  stabilized  at 

each  point.  The  largest  pressure  changes  occur  between  30'  and  80  , which  corresponds 
to  the  changes  in  external  pressures  on  the  APT  as  recorded  during  earlier  flight  tests 
and  noted  in  Figure  3. 

VI . Conclusions  and  Recommendations 


The  primary  objective  of  these  tests  was  to  show  that  the  aerodynamic  environment 
due  to  changing  cabin  pressure  would  not  cause  damage  to  equipment  or  injury  to  per- 
sonnel. Tests  show  that  aerodynamics lly  there  is  no  reason  to  prevent  operation  of 
the  APT  without  output  and  input  windows.  Dynamic  loading  inside  the  APT  is  low,  and 
cabin  pressure  transients  are  relatively  slow  and  will  not  degrade  flight  or  test  crew 
performance . 
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Figure  9.  Dynamic  Pressure  at  Base  of  APT  During  Ground  Test 
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CONVECTIVE  HEAT  TRANSFER  AUGMENTATION  BY 
COMBINED  TECHNIQUES 

"k  kk 

R.  C.  Oliver  and  D.  N.  Black 


Abstract 


Current  Air  Force  developments  in  lasers  and  advanced  propulsive  systems  have  demon- 
strated the  requirements  for  high  rates  of  heat  transfer.  Many  augmentative  techniques 
exist  and  this  paper  reports  on  the  experimental  evaluation  of  a combination  of  two  of 
the  most  promising  convective  techniques,  nucleate  boiling  and  surface  promoters. 

Results  indicate  that  these  techniques  can  be  used  in  concert,  with  additive  effects;  no 
detrimental  effects  were  noted  in  the  range  investigated. 

I . Introduction 

Current  Air  Force  and  industry  research  and  development  efforts  have  illuminated  a 
requirement  for  enhanced  heat  transfer  rates.  Several  specific  areas  of  application  are 
laser  mirrors  and  turbine  blade  cooling.  Increased  heat  transfer  techniques  could  also 
be  employed  in  many  standard  cases  resulting  in  decreased  heat  exchanger  size  or  increased 
performance  or  efficiency.  Successful  developments  of  some  advanced  systems  may  require 
the  use  of  augmentative  heat  transfer  devices. 

The  baseline  case  for  convective  heat  transfer  is  a constant  area,  smooth,  straight 
channel  with  only  gravitational  forces  present.  The  multitude  of  available  techniques 
are  found  in  several  publications,  those  by  Bergles  (Ref.  1 and  2)  being  perhaps  the 
most  complete.  Very  briefly  the  techniques  range  from  passive,  where  some  surface  treat- 
ment is  considered,  to  active,  where  mechanical  devices,  vibration,  injection  or  suct- 
ion and  electro-magnetic  fields  are  employed.  Additionally,  devices  such  as  swirl  flow, 
surface  tension  or  additives  comprise  another  distinct  class  (displaced  promoters). 

Add  to  this  the  application  of  these  techniques  to  either  boiling  or  non-boiling  proc- 
esses and  you  began  to  see  the  magnitude  of  the  possible  combinations.  Of  course, 
combined  techniques  of  any  two  or  more  of  the  basic  techniques  employed  simultaneously 
should  also  be  considered  while  seeking  the  "optimum"  augmentative  device. 

Bergles  summarizes  the  research  to  1973  by  noting  that  passive  techniques  appear 
most  practical  in  the  near  term  and  that  the  only  active  techniques  having  present 
practical  utility  are  mechanical  aids.  He  also  noted  that  work  on  compound  techniques 
should  be  extended,  as  the  few  available  studies  indicate  an  Increased  effectiveness 
(Ref.  2). 
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X I . Theory 


The  system  selected  for  study  was  a displaced  promoter  with  nucleate  boiling.  Types 
of  displaced  promoters  include  disks,  inserts,  ring  inserts,  etc.  with  the  purpose  of 
tripping  the  boundary  layer.  The  repeated  disruption  of  the  boundary  layer  reduces  the 
boundary  layer  thickness  and  enhances  the  heat  transfer. 

Theoretically  this  effect  is  shown  by  considering  the  local  Nusselt  number  developed 
for  laminar  flow  but  also  applicable  to  turbulent  flow  (Ref.  3)1 

N = .332  Pr 3 Re  (11 

u x 


or 


h 

x 


k -332  Pr 
s''  I'x/U 


(21 


Where  N is  the  local  Nusselt  number,  Pr  the  Prandtl  number  and  Re  the  local 
u x 

x 

Reyn.  Ids  number.  Further,  h,  k,  V , x,  and  U are  the  convective  heat  transfer  coeffic- 
ient, thermal  conductivity,  viscosity,  distance,  and  the  velocity,  respectively.  From 
this  we  note  that  as  the  boundary  layer  distance,  x,  is  decreased  toward  zero,  h tends 
towards  infinity,  which  is  certainly  in  the  right  direction. 


The  validity  of  this  trend  has  been  verified  by  a number  of  researchers  listed  by 
Bergles  (Ref.  41,  while  Webb  and  others  have  developed  predictive  techniques  for  re- 
peated rib  or  ring  geometry  using  non-boiling  flows  (Ref.  5 and  61. 

It  was  our  goal  to  demonstrate  that  this  augmentation  could  be  maintained  in  con- 
junction with  nucleate  boiling  enhancement,  to  verify  that  the  effects  are  additive. 
Additionally,  consideration  of  the  effect  of  the  promoters  must  also  include  changes  in 
the  critical  heat  rate.  Yang  reported  a 15-301  reduction  of  sub-cooled  burnout  heat 
fluxes  due  to  large  scale  roughness  (Ref.  71.  Bergles  noted  that  surface  conditions  do 
not  appear  to  significantly  alter  the  boiling  curve  for  reasonably  high  flow  velocities 
(Ref.  81,  while  Durant  noted  an  increase  in  sub-cooled  critical  heat  flux  with  relat- 
ively course  mechanical  roughness  (Ref.  91.  Cambill,  using  swirl  flows,  reported  on 
increase  of  200"!.  in  the  critical  heat  flux  (Ref.  101.  In  summary,  we  must  eventually 
consider  the  effect  of  such  augmenters  on  the  critical  heat  rate;  this  effect  appears 
to  be  presently  undetermined.  The  experimental  plan  reported  here  includes  investi- 
gations of  convective  flow  with  nucleate  boiling  witli  and  without  the  use  of  a displaced 
promoter.  Hie  critical  heat  flux  was  investigated  but  not  quantified.  Primary  emphasis 
was  in  the  boiling  phase  prior  to  the  critical  heat  rate. 
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III.  Experimental  Design  and  Data  Collection 


The  experiment  was  conducted  in  the  Physical  Chemistry  Laboratory  at  the  U.S.  Air 
Force  Academy  from  October  through  December  1977.  The  test  apparatus  consisted  of  an 
aluminum  finned  surface,  11.43  by  11.43  cm  by  1.9  cm  thick.  Uniform  two  dimensional  fins 
1.27  cm  thick  were  used  on  the  fluid  side  in  the  direction  of  flow.  The  fins  maintained 
flow  direction  and  increased  the  heat  transfer  from  a TEMCO  Hot  Plate  Model  HP-1915B 
mounted  under  the  non-finned  side  of  the  device.  Plexiglass  mounted  across  the  device 
(on  the  fins)  formed  a water  tight  apparatus.  A manifold  at  each  end  introduced  and 
removed  water.  Iron-constantan  thermocouples  were  used  to  measure  the  incoming  and 
exit  fluid  bulk  temperatures  as  well  as  the  device  surface  temperature.  Mass  flow  rates 
were  measured  with  a graduated  cylinder. 


A typical  run  would  commence  at  an  established  mass  flow  and  heat  rate  below  the 
onset  iX^nucleate  boiling.  Heat  transfer  parameters  were  noted  in  the  steady  state  and 
then  the  heat  rate  was  increased.  A number  of  points  were  sampled  10  to  15'  F apart 
until  the  crit'i^al  heat  rate  was  approached. 

Several  configurations  were  tested.  The  baseline  condition  was  the  clean  config- 
uration, as  previously  described.  Modifications  using  wires  forming  rectangular  ring 
inserts  were  installed  and  tesHted.  Densities  of  0.3937  and  0.8  inserts  per  centimeter 
(spacing  of  2.54  and  1.27  cm)  weriNexamined.  Data  collection  for  these  cases  was  as 
described  for  the  baseline  case.  ^ . 
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I V . Discussion 


The  experimental  data  were  reduced  and  plotted  as  the  rate  of  heat  transfer  versus 
the  surface-fluid  bulk  temperature  difference.  Figure  1 presents  the  results  of  all 
three  cases,  forced  convection  with  nucleate  boiling,  clean  and  with  both  insert  con- 
figurations . 

As  the  figure  illustrates,  the  critical  heat  rate  was  reached  at  a smaller  temper- 
ature difference  for  the  two  cases  using  wire  inserts  for  boundary  layer  disruption. 

The  critical  heat  rate  for  the  baseline  case  occurred  at  a temperature  difference 
(surface  - fluid  bulk)  of  175  F.  For  the  lower  insert  density  treatment  the  temper- 
ature difference  was  125  F.  The  higher  density  treatment  resulted  in  temperature  diff- 
erence of  100  F.  From  the  data  analyzed,  no  significant  difference  in  the  critical 
heat  flux  was  observed  for  these  treatments.  If  higher  fluxes  were  used  perhaps  some 
differences  would  be  noted.  However,  the  object  of  this  study  was  heat  transfer  below 
the  critical  flux.  Hie  increase  in  insert  density,  then  results  in  enhanced  convection. 
For  a given  flux  the  decrease  in  the  temperature  difference  is  indicative  of  an  in- 
creased convective  heat  transfer  coefficient. 
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Heat  transfer  coefficient  to  the  maximum  temperature  difference  in  each  case  resulted 
in  values  of  23.71,  41.03,  and  71.20  BTUH  per  Sq  Ft-  F.  The  absolute  values  of  these 
coefficients  is  low  because  of  the  limited  mass  flow  rate.  The  flow  had  to  be  maintained 
at  low  rates  to  insure  reaching  the  critical  flux  with  the  fixed  power  heater.  Hie 
possibility  of  different  results  at  higher  flow  rates  must  be  considered  and  experiments 
should  be  performed  to  confirm  this  result.  Such  determination  was  beyond  the  scope  of 
this  effort. 

V.  Conclus ion 


The  data  collected  indicated  that  the  use  of  displaced  promoters,  namely  rectang- 
ular inserts,  increases  heat  transfer  rates  for  flows  where  nucleate  boiling  is  present. 
For  systems  using  nucleate  boiling,  then,  the  additional  use  of  inserts  will  Increase 
the  heat  transfer  rates.  No  mention  of  the  increased  pressure  loss  is  made,  but  its 
effect  must  also  be  considered.  It  is  also  highly  likely  that  even  a third  or  fourth 
technique  could  be  coupled  with  these  to  further  increase  the  heat  transfer  rate. 
Specifically,  fluid  additives  and  surface  treatments  which  increases  the  nucleation 
sites  should  be  considered. 

The  data  collected  and  presented  here  imply  that  further  increases  in  the  insert 
density  will  further  increase  the  heat  transfer  rate.  While  this  is  true,  basei  on  both 
theory  and  the  reported  results,  the  rate  of  increase  decreases  with  increasing  insert 
density  (diminishing  returns).  The  optimum  design  would  have  to  consider  the  increased 
pressure  losses,  initial  cost  and  operating  expenses,  perhaps  more  frequent  defouling. 
The  combined  techniques  do,  however,  have  additive  effects  and  can  be  used  in  concert 
to  further  increase  heat  transfer  rates. 

The  data  collected  are  not  deemed  sufficiently  general  in  flow  and  heat  transfer 
rates  to  justify  empirical  predictive  equations.  The  data  do,  however,  indicate  the 
possibility  of  using  these  two  techniques  in  concert  and  encourages  the  research  into 
the  use  of  additional  techniques  with  these  or  other  combined  techniques. 
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LASER  INDUCED  HEATING : PREDICTION  OF  FRONT  SURFACE  MELTING  TIMES* 

v. 

Raymond  A.  Shu  Is  tad 

Abstract 


A critical  thickness  for  a slab  subjected  to  laser  radiation  is  defined.  It  is  then 
shown  that  for  slabs  of  thickness  greater  than  this  critical  thickness  a simple  equation, 
derived  from  a half-space  solution  to  the  heat  equation,  may  be  used  to  predict  front 
surface  melting  times.  In  addition  it  is  shown  that  for  thin  slabs,  of  thickness  less 
than  three-fifths  of  the  critical  thickness,  another  simple  equation  may  be  used  to 
compute  melting  times.  Further,  front  surface  melting  times  for  slabs  of  intermediate 
thickness  are  shown  to  be  bounded  by  the  times  computed  from  the  two  equations  cited 
above . 


1 . 1 ntroduc  tlon 


An  Important  parameter  in  any  investigation  of  the  effects  of  high-power  laser  radt 
atlon  on  solids  is  the  time  required  to  melt  the  front  surface.  Two  elementary  solut- 
ions of  the  heat  conduction  equation  as  given  by  Garslaw  and  Jaeger  (lief.  11  are  parti- 
cularly useful  for  predicting  surface  thermal  response  up  to  the  time  when  the  front 
surface  begins  to  melt.  1'hese  solutions  assume  one  dimensional  heat  flow  and  neglect 
radiative  and  convective  losses. 


fhe  front  surface  temperature  of  a half-space  of  thermal  conductivity  k and  dlffns- 
lvit\  ■ initially  at  temperature  T(  and  subjected  to  a flux  of  uniform  intensity  1 with 
a fraction  n absorbed  at  the  surface  is: 


surf 
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the  time  to  initiate  front  surface  melting,  7 , is  calculated  in  equation  2: 
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where 
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If  a critical  thickness,  1 , , is  defined  as: 
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then  equation  2 nuiy  he  rewritten  as; 
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It  will  be  shown  In  this  paper  that  If  a slab  Is  of  thickness  creator  than  i , , then 

crlt 

the  time  required  to  melt  its  front  surface  is  also  given  by  equation  4 (i.e.,  the  same 
as  predicted  by  the  half-space  equation!. 

For  uniform  absorbed  flux,  0(1,  on  the  surface  of  a homogeneous  slab  ol  thickness  (, 
the  front  surface  temperature  is  given  by  : 
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Equation  5 is  transcendental,  and  an  iterative  procedure  must  he  us  d if  melt  times  are 
desired.  For  thin  stabs  (<*‘  t/i*  > > 11  equation  5 can  be  reduced  to 

T , * T + (6) 

surf  o k*  3k 

from  which  the  time  required  to  initiate  front  surface  melting,  T.j,  is  given  by 
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Using  equation  3,  Tgj  is  also  found  to 


T r 1 crft  - i (Ss 

Tsl  2.2*:  3*:  <8’ 

It  is  interesting  to  note  that  the  first  term  in  equation  7 is  the  melting  time  for  a 
very  thin  slab  under  uniform  heating  conditions  (i.e.,  dT/dx  » 01. 

Thus,  It  is  seen  that  front  surface  melting  times  are  easily  calculated  for  "thick" 
slabs  via  equation  4 and  for  "thin"  slabs  via  equation  8,  but  one  must  work  considerably 
harder  to  determine  melting  times  for  slabs  of  intermediate  thicknesses.  Ails  investi- 
gation establishes  the  range  of  applicability  of  equation  4 and  8 in  terms  of  the 
critical  thickness,  jl  defined  in  equation  3 and  derived  in  the  next  section  of  this 

paper.  It  will  also  be  shown  that  melting  times  as  predicted  by  equations  4 and  8 may 
be  used  to  bound  the  actual  melting  time  of  slabs  of  Intermediate  thicknesses. 

11.  Definition  of  Critical  Thickness 

An  iterative  procedure  was  developed  to  compute  front  surface  melting  times  for  slabs 


l ?! 


of  various  materials  and  thicknesses  (Fq  Si  and  Intensities.  These  times  were  compared 
to  r , the  time  required  to  melt  the  front  surface  of  a half-space  suhterted  to  the  samp 
flux  U'q  21.  Figure  l depicts  typical  results  obtained  In  these  studies.  The  results 
of  the^e  comparisons  suggested  the  existence  of  a "critical  thickness"  which  Is  defined 
to  he  the  minimum  thickness  wit  ere  the  time  required  lo  melt  the  front  aurtace  ol  the  slsb 
Is  equal  to  the  time  predicted  hy  the  halt-space  solution  (pi. 


With  this  definition  in  hand.  Information  as  to  what  factors  determine  critical  thick 
ness  was  sought.  Hy  requiring  the  front  surface  of  the  slab  to  melt  at  time  equal  r,. 

It  was  found  that 


seemed  to  he  a fundamental  dimensionless  parameter.  Utilisation  ol  critical  thickness 
data  acquired  In  the  comparison  study  described  above  showed  that 


crlt  Or  1 


where  z - 2.2. 


Critical  thicknesses  and  correspond Inn  "*  values"  tor  a varltv  of  net  a l s are  presented  in 
table  1.  These  data  show  that  z Is  Independent  of  material  properties.  Further  studies 
verified  that  this  parameter  Is  also  Independent  of  absorbed  Intensity. 

Thus,  it  Is  proposed  that,  for  slabs  of  thickness  greater  than  the  critical  thick- 
ness (Kq  3),  one  may  with  confidence  predict  front  surface  melt  inn  times  via  equations 
2 or  4.  Further,  It  Is  suggested  that  p may  be  used  as  an  upper  bound  on  front  surface 
melting  times  for  slabs  which  do  not  meet  this  criterion. 

It  was  pointed  out  above  that,  for  "thin"  slabs,  equation  7 or  8 can  be  used  to  pre- 
dict front  surface  melting  times.  In  the  course  of  this  study,  It  was  found  that  "thin 
slab"  could  also  be  defined  In  terms  of  the  critical  thickness.  With  this  definition. 

It  will  be  shown  that  the  range  of  applicability  of  equation  7 or  8 can  be  established. 

111.  Ue  fin  it  lor.  of  Thin  Slabs 

Using  equations  i and  7,  It  la  easy  to  show  that  the  front  surface  melting  tins'  for  a 
slab  ts  given  by 

r.r.i  + f OH 


where  t*  - L'  \ 
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For  thin  slabs  the  aerie.-  converges  extremely  rapidly,  and  t‘  mav  he  adequate'y  approx- 
imated by  the  first  term.  lima. 


A thin  alab,  for  ovir  purpose , amv  be  defined  as  one  thin  enough  such  that  the  series 
contribution  to  the  surface  thermal  response  at  melting  Is  negligible.  This  enables  one 
to  use  the  very  simple  equation  7 to  predict  front  surface  melting  times  and  requires 


where  f Is  some  small  fraction  which  we  arbitrarily  choose  to  be  0.01.  Since  Is  a 

lower  bound  for  T,  equation  13  Is  satisfied  if 


Using  the  expression  for  T ^ j given  previously  (Kq  bf,  equation  IS  may  he  rewritten  as  : 
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It  Is  easily  shown  that  equation  lb  is  satisfied  when 
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thus . It  Is  proposed  that,  for  slabs  of  thickness  less  than  about  three-fifths  of  a 
critical  thickness,  equation  7 or  8 may  be  used  to  predict  front  surface  melting  times 
(l.e.  , t f j),  Further,  since  the  series  contribution  Is  neglected  In  equation  7, 
y may  be  used  as  a lower  bound  on  front  surface  melting  times  for  slabs  of  intermediate 


thicknesses  l.b2  ( 
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IV.  Conclusion 

This  paper  defines  a critical  thickness  for  a slab  subjected  to  a constant  flux  of 
laser  radiation.  this  critical  thickness  can  be  used  In  predicting  front  surface 


molting  time  bounds  for  slabs  of  intermediate  thicknesses.  l'ablt*  2 sutnnari/es  tin-  re- 
sults ot  this  study  In  terms  ol  the  critical  thickness  defined  by  equation  J. 
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Table  1 

CRITICAL  THICKNESS  VS  Z FOR  VARIOUS  METALS 


Metal 

k 

J/cm-sec-k 

At* 

k m 

r ** 

^crlt 

z 

Aluminum 

1.84 

611. 

0.51 

2.27 

Magnesium 

0.96 

605. 

0.26 

2.22 

Titanium 

0.145 

1600. 

0.11 

2.38 

Stainless-Steel 

0.240 

1400. 

0. 16 

2.38 

Lead 

0.320 

300. 

0.04 

2.08 

Copper 

3.91 

1050. 

1.50 

1.82 

Nickel 

0.67 

1428. 

0.46 

2.38 

Cadmium 

0.944 

294. 

0.12 

2.17 

Silver 

4.18 

935. 

1.50 

1.92 

r • 2.18  + .06 

- 300  0 k 

> 

1 - 5 x 10J  w/cm* 


Table  2 


PREDICTION  OF  FRONT  SURFACE  MELTING  TIMES 
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'.ill'  PHI  I'kl  l'lCAI.  BASIS  ANT>  l WU  MKNTAl.  APPLICATION 
or  VUE  ANALOGY  RKtVFKX  MASS  AND  HEAT  TRANSFER 

S.  1 . I card! 


Abstrac  t 

This  paper  outlines  the  theoretical  basis  for  the  analogy  between  mass  and  heat  trans- 
fer. It  discusses  an  experimental  method  for  estimating  temperatures  in  flows  by  relating 
them  to  the  measurement  of  mass  concentrations  of  a trace  gas  in  the  flow. 

I.  Introduc  t Ion 

Due  to  concern  over  the  impingement  of  hot  exhaust  gases  from  the  gas  dynamic  laser 
on  the  fuselage  and  control  surfaces  of  the  Airborne  laser  Laboratory.  an  experimental 
technique  was  developed  to  map  the  exhaust  plume  and  estimate  its  temperature.  Vhe 
experiments  Involved  modeling  the  flow  in  the  wind  tunnel.  Introducing  a trace  gas  Into 
a cold  let.  and  then.  measuring  the  mass  concentrations  of  this  trace  gas  at  various 
locations  on  the  surface  of  the  aircraft  model.  By  applying  the  analogy  between  mass 
and  beat  transfer  an  estimate  of  the  exhaust  plume  temperature  could  be  made.  litis  paper 
outlines  the  theoretical  basis  for  these  experiments. 

II.  >e  fin  it  ion  of  Verms 


Lewis  Number  he. 


- dvnamtc  viscosity  coefficient 

Pr  * brand t l Number  » Ho- 


diernal diffusion  coefficient 


specific  heat  at  constant  pressure 
thermal  conductivity 


Sc.  Schmidt  Number 


mass  diffusion  coefficient  ,0  D. 


P - density,  mass  unit  volume 

I)  mass  diffusion  coefficient  for  component 

■*  j in  a multi-component  mixture 

a molecular  thermal  diffusivlty  k/oc 
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of  heat  equal  to  rocT  and  fluid  flowing  downward  carries  an  amount  m'  cl"  if  I"  since  m 
in'  the  heat  flux  per  unit  area  through  aa  is: 

me  if-T  ) , . 
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The  turbulent  exchange  causes  a mass  flow  rate  per  unit  area  of  one  component  of  the 
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Dividing  heat  flux  by  mass  flow  race  we  have  the  following  for  turbulent  flow 


For  the  analogy  to  hold  between  laminar  and  turbulent  flows  then 


a JL. 
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We  are  now  able  to  relate  the  mass  fraction  difference  to  the  temperature  difference 
in  a flow.  If  the  Lewis  number  is  one,  we  have  seen  above  that  the  ratio  of  heat  flux 
to  mass  flux  at  any  particular  point  in  the  flow  is  the  same  regardless  of  the  type  of 
flow  (laminar  or  turbulent).  It  follows  that  the  relationship  between  temperature  and 
mass  concentrations  at  any  particular  point  in  the  free  stream  and  any  two  other 
locations  in  the  flow  are  also  the  same. 


adjusting  this  equation. 
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Experimentally , If  we  desire  to  estimate  the  temperature  In  a flow  (knowing  the 
temperature  In  the  free  stream  and  at  one  other  location  In  the  flow),  we  can  do  so  by 
modeling  the  flow  and  Introducing  a trace  gas  Into  the  flow.  Measuring  the  mass  fraction 
of  the  trace  gas  at  these  locations  and  using  the  following  equation,  we  can  estimate  the 


temperature 


(T  - T ) + T 


This  equation  of  course  assumes  that  the  mass  fraction  of  the  trace  gas  in  the  free 
stream  Is  zero  (01^  » 0) 

IV.  Conclvjlon 

Experiments  verifying  this  technique  have  been  performed  by  a number  of  researchers 
and  are  sighted  In  Reference  2 pg.  472.  Tills  technique  was  used  to  estimate  temperatures 
In  hot  exhaust  gases  of  the  Airborne  l,aser  laboratory  I.aser  Exhaust  Plume.  Tire  experi- 
ments were  run  In  the  subsonic  wind  tunnel  at  the  U.S.  Air  Force  Academy.  The  results 
of  these  experiments  are  reported  in  Reference  3. 
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Ahat  l ap  1 


llils  papor  dosprlhos  an  alt  ornat  tvo  olaotronlc  1111  and  drag  measuring  svatpm  for 
tlio  Popartmoiit  ol  Aeronautics  l«’"  wind  timin' Is.  An  analvsta  ol  flip  problem  areas  of  tlio 
present  median  Ira  l measuring  svstom  Is  given.  An  alternate  it  1 1 1 m 1 oloetronlr  svstom  Is 
described  along  with  flip  lost  vpsnlts  obtained  In  tost  tng  this  a 1 tenia  ft ve,  Kecommend - 
at  Ions  ol  tntnvp  changes  sml  1 mp  1 i-nn'ii  t • t Ion  ot  this  now  svstom  aro  made. 


I.  lnt  voilnp  t li'ii 

t'ommerclal  Iv  available  wind  t Illinois,  with  a 1.'"  ill  amp  tor  pironlar  tost  soil  Ion,  aro 
napil  hv  tin-  l'p\>artiiipiit  ol  Aorouant  ios  at  flip  Ihilti'il  Statos  All  t'orco  Aoademv  for  demon  - 
st  rations  anil  experiments  In  Ini  roiluot orv  Aeronautical  tnglnoortng  ponrsos  . lliov  aro 
mpi- haul  pal  Iv  slin|'lp,  ami  falvlv  tnoxpoiis  1 vo . hut  aro  papal'  lo  ol  producing  ilaia  lot 
ipiant  It  at  l vs'  ana  l vs  os . 

Hip  wlml  tiiniiola  aro  oonstrnotoil  with  a mpphanlpal  pointer  ami  apato  svstom  to  ills 
plav  lilt  ami  drag  lorpos  apt  lug  on  tlio  tost  moilol.  Usual  Iv,  howovor,  tlio  pointers 
lliiptuatp  excessively,  wlilph  makes  data  polloptlon  illfllpiilt  and  iinportain.  Ilils  papor 
reports  on  tin-  feasibility  ,>l  Improving  flip  present  mpphanlpal  measurement  svstom  nsod 
In  p lassrootn  wind  t mmols  hv  Installation  ol  a now  oloptronlo  dlpltal  inoastiromont  svstom. 

11  io  design  ol  tlio  linprovod  a It  prnat  1 ve  prosontod  lioro  was  groat  Iv  influenced  hv  tlio 
criteria  that  tlio  svstom  ho  t nospons  1 vo  and  oasilv  tnstallpd.  Addtttonallv  all  olop- 
tronlp  plrpultrv  required  had  tv'  ho  ponst  mot  od  hv  an  undergraduate  stmlont  using  roadllv 
aval  I ship  parts. 

II.  Wind  I'nnnol  t’onst rno t ton 

figure  I shiiws  tho  ponst  nipt  ton  ol  tlio  olaasroom  wind  tinmols  nsod  at  tlio  PSAV 
Academy,  llio  wind  tinmols  havo  a l.'  Inph  diameter  tost  soot  ton.  A model  ts  plaood  on  a 
tost  stand  wlilph  Is  oonnootod  to  inoohaiitoal  potntors  that  tndlpato  tho  amount  s of  lilt 
ami  drag.  Air  ontora  tlio  tmmol  through  tho  honovpomh  Inlot.  llio  honovpomh  forms  a 
sortos  ol  small  tnhos  through  wlilph  flip  air  Is  drawn,  straightening  and  smoothing  tho 
airflow.  tin'  wall  ot  tlio  tost  soot  ton  ts  madp  ol  oloai  Plexiglass,  onahling  tho  mods  l 
to  ho  vlpwod  dining  tosting.  Anglo  ot  at  lark  nil  pan  ho  adjustod  hv  a t Immh  -screw  on 
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the  test  stand  through  th<*  range  -ill  to  +20  . 
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K 1 gore  1.  Wind  Tonne 1 Construct  loti 

Zeroing  thisnb-screws  ar*>  provided  to  set  the  pointers  to  zero  under  no-wlnd  conditions. 
Air  Is  drawn  through  the  tunnel  by  a four-bladed  electric  lan  mounted  at  the  wind  tunnel 
exit. 

111.  Mechanical  Korea-Measuring  System 
A.  Lift  and  brag  vs  Angle  of  Attack 

For  low  auglea-of -attack , C(  change*  nearly  linearly  with  angle  of  attack,  while 
the  c curve  Is  parabolic,  brag  la  alwavs  positive,  whether  a Is  positive  or  negative, 
while  lilt  can  take  on  positive  or  negative  values.  Figure  2 shows  an  ideal  plot  of 
lift  and  drag  coefficients  vs  Or. 
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Figure  2.  Typical  I lit  and  Drag  Coefficients  vs  a 
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1'.  Pointer  Construction 
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Figure  3 shows  the  present  method  of  measuring  and  Indicating  lift  and  drag 
forces  with  the  classroom  wind  tunnels.  Lift  and  drag  are  measured  using  the  same  method 
(only  the  1 1 ft -measuring  system  Is  shown 1 . Lift  forces  the  model,  test  stand,  and  steel 
spring  upward,  pulling  a string  In  the  direction  Indicated.  The  string  pulls  on  the 
pointer,  which  is  balanced  by  a counterweight  and  a small  spring.  llie  small  spring  Is 
omitted  from  Figure  ' for  simplicity.  As  lift  Increases,  the  string  Is  pulled  and  the 
pointer  moves  toward  the  right.  Indicating  an  Increase  in  force. 


Mgure  3.  Mechanical  Force  Measuring  System 

As  lift  decreases,  tension  In  the  string  Is  relaxed,  and  the  counterweight  and  steel 
spring  force  the  pointer  hack  to  the  left,  Indicating  a decrease  In  force.  The 
mechanical  system  Is  not  capable  of  Indicating  downward-acting  lift. 

I'.  Pause  of  Pointer  Fluctuation 


Phe  steel  spring  Is  not  very  stiff,  and  the  pointer  assembly  Is  of  low-mass 
construction.  Turbulent  flow  Inside  the  wind  tunnel  causes  the  model  to  vibrate  con- 
siderably. The  fan  causes  the  entire  wind  tunnel  to  vibrate.  Together,  these  two 
effects  cause  the  pointer  to  fluctuate  erratically,  making  accurate  readings  very 
difficult. 


D.  Resultant  nata  and  Uncertainty 

Most  of  the  test  models  used  In  this  particular  wind  tunnel  (1/72  acale  models') 
produce  a maximum  ltfttng  force  of  0.70  pound.  Some  larger  models  (1/4F  scale)  produce 
up  to  one  pound  of  lift.  Typical  values  of  lift  change  by  only  0.02  pound  when  Of  la 
changed  by  two  degrees.  Corresponding  changes  In  values  of  drag  are  only  about  0.01 
pound.  However,  the  pointers  fluctuate  up  to  0.06  pound  for  both  lift  and  drag.  This 
fluctuation  covers  a greater  range  than  the  distance  between  two  data  points.  A 0,06 
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lift  COEFFICIENT  (CL) 


pound  fluctuation  for  a model  producing  0.70  pound  of  lift  equates  to  an  uncertainty  of 
up  to  97..  This  same  fluctuation  for  a maximum  drag  of  0.22  pound  (typical  for  a 1/72 
scale  model)  equates  to  an  uncertainty  of  up  to  37%.  Figure  4 shows  a typical  range 
of  pointer  fluctuation  against  a calibrated  scale. 


Seale 

Mange  of  Fluctuation 

Pointer 


Figure  •*.  Pointer  Fluctuation 

When  data  Is  taken  and  lift  and  drag  coefficients  are  plotted  against  angle-of- 
attack,  a large  range  of  uncertainty  Is  evident  as  Is  shown  In  Figure  5.  It  is  desir- 
able then  to  reduce  the  range  of  uncertainty  as  much  as  possible,  l'he  Improved  Instru- 
mentation system  presented  here  was  developed  to  accomplish  this  goal. 


Figure  S.  Mechanical  Data  Uncertainty 
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I V . Digital  Measurement  System 


l'he  use  of  electronic  circuitry  to  solve  the  pointer  fluctuation  problem  was  chosen 
because  of  Its  simplicity  and  versatility.  Voltage  fluctuations  caused  by  the  wind 
tunnel  vibrations  can  easily  be  smoothed  into  steady  values  with  filters,  and  electronic 
signals  can  be  easily  controlled  and  measured. 


A.  Force  Sensing 


Strain  guagcs  were  used  to  sense  lift  and  drag  forces.  These  devices  measure 
the  amount  of  bending  that  the  steel  springs  experience.  An  enlarged  strain  gauge  is 
shown  In  Figure  b.  The  foil  pattern  has  a nominal  resistance  of  120  ohms.  However, 
when  the  strain  gauge  Is  stretched.  Its  resistance  will  Increase  by  a slight  amount. 


Solder  Point  - 
Resistive  Foil 


Plastic  Tacking 


Mounting  Guide 


Conversely,  the  resistance  will  decrease  slightly  if  the  strain  gauge  Is  compressed. 
The  change  in  resistance  Is  proportional  to  the  amount  of  stretch  or  compression  the 
strain  gauge  experiences.  The  strain  gauges  were  mounted  on  the  surfaces  of  the  steel 
springs  as  shown  In  Figure  7. 
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Figure  7.  Strain  Gauge  Mounting 


In  this  configuration,  as  a spring  Is  deflected  downward,  the  strain  gauges  would  be 
stretched,  and  their  resistances  Increased.  On  the  other  hand.  If  the  spring  were  de- 
flected upward,  the  strain  gauges  would  be  compressed,  and  their  resistances  decreased. 
Figure  8 Illustrates  the  stretching  and  compressing  of  the  strain  gauges  upon  bending  of 
the  springs. 
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Strain  Gauge  Stretch  and  Compression 


B.  Strain  Gauge  Mounting  Problems 


Strain  gauges  must  be  precisely  mounted  to  give  proper  Indications  of  spring 
deflection.  Normally,  a reference  line  is  scribed  on  the  surface  to  which  the  strain 

and  the  gauge  is  aligned  with  this  reference  line,  as  shown 
However,  because  the  test  stand  springs  were  located  in  hard-to-reach 
locations,  a reference  line  could  not  be  scribed. 


gauge  is  to  be  mounted 
in  Figure  9 
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ideal  Strain  Gauge  Mounting 


F igure  9 


As  a result 


As  near  proper  as  possible  strain  gauge  alignment  was  estimated  by  eye 
the  strain  gauges  ended  up  slightly  misaligned  as  shown  in  Figure  10. 
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Figure  10 


To  properly  mount  the  strain  gauges  In  any  future  endeavors,  the  spring  unit  should 
be  removed  from  the  wind  tunnel.  If  this  were  done  the  springs  could  then  be  placed  on 
a suitable  workbench  where  proper  mounting  procedures  could  be  followed.  Tills  was  not 
done  In  this  first  modification  because  of  the  time  Involved  and  the  difficulty  of 
disassembling  the  wind  tunnels. 

C.  Conversion  of  Force  to  Electrical  Signal 

Through  the  use  of  these  strain  gauges,  lift  and  drag  forces  on  the  model  pro- 
duced an  electrical  resistance  In  an  electrical  circuit.  The  strain  gauge  resistances 
changed  linearly  with  the  amount  of  force  acting  on  the  model.  The  gauges  were  arranged 
In  a Wheatstone  Bridge  circuit  as  shown  in  Figure  11.  A change  in  aerodynamic  force 
caused  a proportional  change  in  voltage.  The  potentiometer  was  adjusted  to  provide 
zero  output  voltage  under  no-force  conditions. 


Voltare 

Source 


Figure  11.  Strain  Gauge  Bridge  Circuit 

D.  Temperature  Compensation 

Only  two  strain  gauges  were  mounted  to  each  steel  spring  in  the  wind  tunnel. 

The  other  two  strain  gauges  in  each  bridge  circuit  were  dumny  gauges,  mounted  on  a strip 
of  aluminum  and  placed  In  close  proximity  to  the  active  gauges.  As  a result,  all  four 
gauges  were  always  subjected  to  the  same  ambient  temperatures. 

(In  many  circuits  using  active  and  dumny  strain  gauges,  the  active  gauges  are 
remotely  mounted,  while  the  dumny  gauges  are  mounted  with  other  circuitry.  The  two 
types  of  gauges  are  then  subjected  to  unlike  temperatures,  and  they  have  different 
nominal  resistances.  As  a result,  the  electrical  characteristics  of  the  circuit  change 
aa  temperature  changes.  To  avoid  this  problem,  the  dummy  strain  gauges  were  co-located 
with  the  active  strain  gauges.) 

E.  Electronic  Filtering  and  Calibration 

The  strain  gauge  output  signal  with  a 5 volt  power  supply  was  about  1 milli- 
volt (»V)  per  pound  of  force  acting  on  the  model.  However,  the  signal  fluctuated  as 


the  springs  vibrated.  The  signal  had  to  be  filtered  and  amplified  before  the  true 
signal  could  be  accurately  measured.  Figure  12  shows  the  voltage  from  the  strain  gauges 
with  these  random  fluctuations  over  a period  of  time. 


voltage 
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Figure  12.  Unflltered  Strain  Gauge  Voltage 

This  signal  could  not  be  displayed  as  a fixed  value.  However,  the  fluctuations  were 
easily  filtered  with  the  proper  electronic  circuitry.  The  voltage  then  appeared  smooth 
and  took  on  the  average  value,  as  shown  In  Figure  13.  Hie  voltage  was  also  amplified 
to  a level  which  corresponded  to  the  amount  of  force  acting  on  the  model,  so  that  for 
example,  one  pound  of  force  produced  one  volt,  1.5  pounds  produced  1.5  volts,  etc. 


voltage 


Figure  13.  Filtered  Strain  Gauge  Voltage 

A low  cost  operational  amplifier  (op-ampl  using  a feedback  resistor  and  a feedback 
capacitor  was  used  to  accomplish  the  filtering  and  amplification.  The  op-amp  circuit 
is  shown  In  Figure  14.  Amplification,  or  gain,  was  adjusted  by  varying  the  value  of 
feedback  resistor  R_. 
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Figure  14.  Op-Amp  Circuit 
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Gain  for  the  circuit  in  Figure  14  is  given  by  the  equation:  Circuit  gain  * — - — • Filt- 

Ks 

ering  was  accomplished  by  the  parallel  combination  of  and  Cp.  The  time  constant 
formed  by  this  circuit  is  the  product  of  Rj,  (ohms')  X C„  (farads')  and  is  expressed  in 

seconds.  A larger  time  constant  will  allow  the  circuit  to  pass  only  lower  frequencies. 
Fherefore,  it  was  desirable  to  have  a very  large  time  constant  so  that  only  direct 
current  (DC)  was  passed.  Consideration  for  values  of  R^  and  Cp  are  discussed  in 
Section  111.  J. 

F.  Digital  Force  Display 

To  read  the  voltage  produced  by  the  op-amp,  a digital  voltmeter  (DVM)  was  used. 
If  the  amplifier  is  calibrated  to  produce  one  volt  per  pound  of  force,  and  there  are 
1.34  pounds  of  force  acting  on  the  model,  the  digital  display  would  indicate  1,3^0 . 

G.  Lift  Measure  at  Negative  Angles -of -At tack 

rhe  mechanical  pointers,  as  they  come  with  the  wind  tunnels,  are  rigged  so  chat 
they  can  only  indicate  lift  in  an  upward  direction,  so  the  model  can  only  be  tested 
for  lift  at  positive  values  of  «.  However,  the  test  stand  is  adjustable  trom  20  nose- 
up  to  20  nose-down.  Fortunately,  the  digital  measurement  system  produced  a negative 
voltage  which  was  proportional  to  lift  acting  in  the  downward  direction,  and  the  DVM 
used  had  the  capability  to  indicate  negative  voltages;  therefore,  an  added  capability 
for  negative  lift  measurements  was  achieved.  For  example,  if  the  test  model  was 
oriented  in  the  relative  wind  at  or  * -10  , giving  rise  to  a 0.33  pound  of  force  in  the 
downward  direction,  the  Dill  would  display 


H.  Digital  Measurement  Uncertainty 

The  digital  instrumentation  circuitry  was  calibrated  to  produce  10  millivolts 
per  pound  of  force.  The  amplification  was  kept  low  to  reduce  the  cost:  of  the  circuitry; 
die  DVM  used,  however,  was  precise  when  measuring  low  voltage  signals.  The  DVM  display 
fluctuated  within  a range  of  0.010  pound.  For  example,  u 0.43  pound  of  force  would 
cause  the  DVM  display  to  fluctuate  from  "0.425"  to  "0.435"  over  a period  of  time.  This 
uncertainty  is  six  times  less  than  the  0.0C  pound  uncertainty  experienced  with  the 
mechanical  system.  Die  digital  system  produced,  therefore,  1/b  the  data  point  scatter 
of  the  mechanical  system.  Figure  15  on  the  following  page  shows  graphically  how  the 
digital  system  reduced  data  point  scatter.  Digital  uncertainty  at  maximum  values  of 
lift  and  drag  was  under  1.5",  for  lift  and  under  5.0";.  for  drag. 
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Figure  15.  Digital  Data  Uneertalntv 


1.  Digital  Circuit  equipment  Description 
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rhe  external  features  of  the  Filter 'Amplifier  are  shown  in  Figure  16. 


Figure  It.  Filter 'Ampltt ler 

Refer  to  the  Appendix  A for  the  circuit  schematic.  Five  volts  were  applied  across  the 
strain  gauge  bridge  circuits'.  l'o  zero  under  no-force  conditions  5k  f)  potentiometers 
were  used  to  set  the  bridge  circuit  output  voltages.  Zero,  g had  to  be  accomplished 
because  models  of  different  weight  deflected  the  steel  springs  by  different  amounts 
when  they  are  mounted  to  the  test  stand.  Hie  potentiometers  were  10-turn  types  to 
make  zeroing  easier,  as  this  was  a sensitive  circuit. 

rhe  strain  gauge  bridge  circuit  acted  as  a variable  voltage  divider.  As  the  act- 
ive gauges  were  deflected,  the  bridge  circuit  output  voltage  changed.  The  selector 
switch  selected  lift  or  drag  strain  gauge  voltages.  This  voltage  was  then  applied  to 
the  op-amp  inputs.  The  positive  strain  gauge  signal  was  applied  to  the  Invert  input, 
so  the  output  was  180  out  of  phase  with  the  input.  For  positive  deflections,  the 
lift  gauges  were  in  compression,  and  the  drag  gauges  were  in  tension.  For  this  reason, 
duniny  and  active  strain  gauge  arrangements  were  opposite  in  the  lift  and  drag  circuits. 
The  calibration  resistors  (500k  0 multiturn  trim  potentiometers'!  were  selected  by  the 
selector  switch  in  conjunction  with  the  bridge  circuit  voltage  selection.  Varying  the 
values  of  these  resistors  changed  the  gain  of  the  amplifier.  Therefore,  amplification 
could  be  calibrated  differently  for  lift  and  drag  using  only  one  op-amp.  This  kept 
parts  costs  to  a miulnnsn.  Tl»e  1 fi  F capacitor,  in  parallel  with  the  500k  0 resistors, 
formed  an  RC  active  filter  with  the  op-amp.  The  DC  Offset  potentiometer  was  a multi- 
turn  trim  pot.  It  was  adjusted  so  that  the  op-amp  produced  no  output  voltage  when  no 
input  voltage  was  present. 
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J.  Filtering  vs  Zeroing  Problems 


Theoretically , the  digital  system  should  be  capable  of  filtering  all  of  the 
strain  gauge  voltage  fluctuations  as  shown  In  Figure  12.  However,  because  the  airspeed 
fluctuated  and  the  model  and  the  wind  tunnel  vibrated  excessively,  an  excessively  large 
feedback  capcitor  \excessive  In  the  sense  discussed  below)  would  be  needed  to  completely 

smooth  these  fluctuations. 


Originally,  a 5 (iF  feedback  capacitor  was  used,  as  Is  shown  In  Figure  17.  For 
a calibration  of  lOtnV  per  pound,  the  value  of  Rj.  was  set  to  100k  fl. 


500k  * 


Figure  17.  Feedback  Circuit 

l'he  time  constant  was  100k  Ox  5 gV  » 0.5  seconds.  A typical  time  constant  is  shown 
in  Figure  IS.  An  impulse  voltage  at  t > 0 would  cause  the  output  waveform  shown. 


vol  tage 


The  time  constant  produces  a lag  in  response  time  in  the  zeroing  circuit. 

After  the  zeroing  pots  are  stopped  rotating,  the  output  voltage  would  continue  to  rise 
or  fall  for  the  duration  of  the  time  constant.  Oils  time  lag,  coupled  with  the  high 
sensitivity  of  the  circuit,  made  zeroing  a very  difficult  procedure. 

Vo  alleviate  the  situation,  a tradeoff  was  made  between  filtering  qualities  and 
ease  of  zeroing.  A l flF  capacitor  was  substituted  for  the  5 4F  capacitor.  The 
resultant  time  constant  was  100k  Qx  1 UV  * 0.1.  Zeroing  was  made  easier  because  of 


tin 


Che  smaller  time  lug.  There  was  also  very  little  difference  in  filtering  quality 
between  the  two  values  of  capacitors  since  most  of  the  fluctuations  were  of  shorter 
period  than  0.1  second.  (See  Section  V for  an  improved  zeroing  feature.) 

V.  Results 

Figure  19  compares  the  uncertainty  ranges  of  the  mechanical  pointers  and  digital 
display.  The  improved  characteristics  of  the  digital  measuring  system  are  quite 
evident.  The  test  data  showed  that  the  digital  measuring  system  did  reduce  the  un- 
certainty error  by  a large  amount  (6007..  reduction  in  lift  and  as  much  us  a 500/!,  re- 
duction in  drag)  resulting  in  a significantly  improved  dutu  quality. 

VI.  Cost  of  Modification 

The  circuitry  used  in  the  modification  was  designed  to  be  low  cost,  'Die  price  of 
all  circuitry,  excluding  the  DVM,  was  about  $230.  because  the  circuit  used  very  few 
components,  it  had  low  gain.  Output  was  only  about  7mV  for  a 1/72  scale  model  at  max- 
imum lift.  To  read  this  low  voltage  precisely,  an  expensive  DVM  (costing  around 
$1400)  needed  to  be  used.  The  DVM  used  in  this  modification  was  used  by  the  Department 
of  Aeronautics  for  u variety  of  applications  and  therefore  contained  many  features  not 
required  for  this  wind  tunnel  modification.  The  same  digital  accuracy  can  be  main- 
tained by  further  amplifying  the  struin  gauge  output  voltage  from  our  lOmV/lb  to  lV/lb. 
This  higher  voltage  could  then  be  measured  with  a low  cost  DVM,  available  for  about 
$150. 

One  drawback  of  the  initial  circuit  was  that  the  time  constant  changed  when  K^.  was 
adjusted.  This  circuit  was  greatly  Improved  by  using  only  $3.00  worth  of  additional 
parts.  The  improved  circuit  is  shown  in  Appendix  B.  During  zeroing,  the  filter 
capacitor  was  cut  out  of  the  circuit.  This  resulted  in  no  time  lag  during  the  zeroing 
operation,  yet  it  permited  the  cupacltor  to  be  large  enough  for  improved  filtering 
qualities.  The  gain  of  tills  circuit  could  be  set  to  lV/lb.  or  to  lV/Newton.  This 
higher  voltage  could  easily  be  read  to  the  thousandths  place  by  a low  cost  DVM. 

Using  the  improved  circuitry  and  a low  cost  DVM,  the  cost  of  a single  modification 
would  be  $383.  Tile  classroom  wind  tunnels  are  currently  available  at  a cost  of  about 
$2000.  Therefore,  the  cost  of  the  modification  is  only  197.  of  the  total  cost  of  a new 
wind  tunnel. 

VI I . Cone lus ions 

The  classroom  wind  tunnels  In  their  original  configuration  can  only  be  used  for 
basic  aeronautical  demonstrations  with  the  mechanical  measuring  system.  However,  with  the 
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improved  accuracy  of  Che  digital  measuring  system  described  here,  the  wind  tunnels  could 
be  used  for  more  exact  experimentation  by  undergraduate  students. 

The  electronic  circuitry  used  in  the  digltul  measurement  system  was  inexpensive  and 
readily  available  from  electronic  supply  houses.  A student  could  easily  obtain  the 
required  parts  and  build  a new  measurement  system  in  about  one  month. 
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VI XI.  Kaconanenda  tlona 

Ifhe  digital  measurement  system  should  he  reeons  true  Cad  using  t l\o  Improved  circuitry 
provided  (u  Appendix  11.  litis  Improved  circuit  permits  the  use  of  a low  coat  DVN,  If 
there  are  any  future  plana  for  using  the  1.’"  wind  tunnels  for  precise  experiments,  1 
recommend  that  all  wind  tunnels  committed  for  this  purpose  be  modi  fled  to  the  Digital 
display  system. 
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A COMPUTERIZED  UNDERGRADUATE  AIRCRAFT  DESIGN  COURSE 

Hr 

G.  T.  Matsuyama 

ABSTRACT 

The  integration  of  computer  programs  with  an  undergraduate  aircraft  design  course  is 
described  with  emphasis  placed  on  the  computer  program  known  as  AER0464RPV,  which  was 
written  specifically  for  the  course.  This  approach  to  design  education  has  proven  highly 
effective  at  the  Air  Force  Academy  during  the  period  1975  through  1977. 

The  course  is  structured  to  require  periodic  tasks  of  students  throughout  the  semester 
to  Include  (1)  an  initial  weight  estimation,  (2)  wing,  fuselage,  and  tail  design,  and 
configuration  layout,  (3)  lift,  zero-lift  drag,  and  drag-due-to-lif t determination,  (4)  a 
refined  fuel  and  empty  weight  analysis,  and,  finally,  (5)  a final  report.  Manual  calcu- 
lations are  required  for  all  tasks,  and  extensive  familiarity  with  the  design  process  is 
emphasized.  Computer  programs  may  be  used  only  after  initial  manual  calculations  have 
been  completed. 

Two  computer  programs  are  used.  The  program  identified  as  AER0464RPV  was  written 
specifically  for  the  design  course  and  employs  course  methods.  Tills  program  is  discussed 
to  some  extent  in  the  paper.  It  requires  basic  vehicle  geometry,  aerodynamic  coeffic- 
ients, and  a mission  profile  to  be  entered  as  inputs.  Output  is  controlled  by  an  input 
variable,  and  progresses  in  complexity  as  further  Information  becomes  available.  The 
program  computes  weight/ fuel  fractions  and  "flies"  the  vehicle  while  attempting  to 
achieve  convergence  with  a historical  empty  weight.  Early  in  the  semester,  initial 
geometric  data  such  as  wing  characteristics  and  fuselage  volume  are  available  as  output. 
Centers  of  gravity  follow  after  component  locations  are  established.  Fuel  tank  sizing 
is  automatic  and  a refined  weight  analysis  is  possible.  Finally,  fuel  sequencing,  static 
margins,  trim  drags,  and  a time  history  of  flight  are  all  obtainable  by  the  end  of  the 
semester  when  sufficient  vehicle  data  is  available.  This  program's  value  lies  in  the 
fact  that  (1)  only  course  methods  are  employed,  (2)  its  input  and  output  structure  allows 
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Che  student  to  develop  an  Initial  data  deck  and  add  to  It  as  further  Information  Is 
gained,  (31  rapid  evaluation  of  design  changes  Is  possible,  (41  previous  output  may  be 
re-generated  merely  by  changing  one  Input  variable  and  without  having  to  disassemble  or 
rearrange  the  data  deck,  and  (41  Its  structure  allows  for  modification  to  accomodate  any 
fighter-type  vehicle.  The  baseline  vehicle  Is  an  advanced  RPV  using  a .1-85  turbojet 
engine.  To  date,  the  program  has  been  successfully  modified  and  used  to  Investigate  a 
manned,  f Ighter/ tralner/ECM  vehicle,  and  now  Includes  a TFE  731-.'  turbofan. 

The  second  program,  known  as  DIGITAL  DATCOM,  predicts  aerodynamic  coefficients  and 
stability  derivatives  for  the  wing,  tails,  and  fuselage  as  well  as  for  the  Integrated 
configuration.  Its  unique  Input  and  output  capability  allows  the  student  to  Input  either 
one  component  (e.g.  wlngl  and  receive  all  pertinent  Information  or  to  Input  the  entire 
configuration  and  still  receive  Individual  component  results  together  with  the  total 
configuration  results.  The  coefficients  obtained  here  are  used  as  Inputs  to  AKR0464RPV. 
It  should  be  noted  at  this  time  that  DIGITAL  DATCOM  Is  a 1973  version  which  was  heavily 
modified  for  Academy  use.  A new  and  expanded  version  Identified  as  AKFDL-TR-76-45  and 
documented  In  3 volumes  was  recently  published  by  the  Flight  Dynamics  Laboratory,  Wrlght- 
Patterson  AFB,  Ohio. 

To  date,  AER0464RPV  has: 

1.  Successfully  been  employed  over  3 semesters  to  assist  students  In  accompl lshtng 
rapid  design  changes,  and  to  evaluate  the  effects  of  these  changes. 

2.  Demonstrated  a clear  advantage  over  other  computerized  methods  due  to  the  fact 
that  only  course  methods  are  programmed,  thus  allowing  tnnedlate  duplication  of  manually 
acquired  data  and  Instilling  confidence  In  the  student. 

3.  Enhanced  the  student's  learning  process  by  providing  adequate  data  upon  which  to 
base  undergraduate  design  decisions. 

4.  Added  to  the  student's  appreciation  of  the  Interrelationships  between  configur- 
ation, center  of  gravity,  and  aerodynamic  characteristics  due  to  the  bulldlng-up  process 
of  the  data  bank  as  well  as  the  extensive  capabilities  of  the  program. 


5.  Given  Che  student  an  appreciation  of  the  usefulness  of  the  computer  as  a tool 


6.  Demonstrated  an  adaptability  for  conversion  as  evidenced  by  Its  successful  modi 


flcatlon  and  use  where  a multi-role,  manned  fighter  was  investigated 


In  light  of  the  above,  It  is  concluded  that  the  use  of  a baseline  computer  program 


employing  only  course  methods  offers  significant  advantages  to  an  undergraduate  aircraft 


design  course 
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A TECHNIQUE  FOR  VORTICLTY  MEASUREMENT  IN  UNSTEADY  FLOW 

John  K.  Kcesee 
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Michael  S.  Francis 

*** 

.lames  l).  Lang 

AbSTRACT 

A technique  Is  described  for  determining  spatial  vorticlty  distributions  In  non- 
statlonarv  fluid  flows  exhibiting  a periodic,  time-varying  mean  velocity  field.  Although 
vorticlty  Is  one  of  the  roost  significant  physical  variables  In  many  complex  flows  involv- 
ing shear.  It  is  a quantity  which  has  been  difficult  to  measure  accurately.  The  method 
discussed  In  the  present  paper  is  one  based  on  a simple  premise  which  involves  no 
assumptions  on  the  dynamic  nature  of  the  flow  and  is  therefore  inherently  accurate. 

The  basis  for  the  method  lies  in  the  generalized  definition  of  circulation  as 
fol lows : 


T(t)  tf(x,t)  • d7(x) 


r is  the  circulation  about  a closed  spatial  contour,  c. 

U • d*  is  the  inner  product  of  the  velocity  vector  with  „the  local  patli  vector 
on  the  contour,  c. 

Equation  (11  can  be  related  to  the  vorticlty  variable  through  Stoke' s theorem: 


£ U( x\ t ) • d«*(x)  = JJ  • d t(x) 


£ * £(x,y,z,tl  is  the  vorticlty  vector. 

A is  a capping  surface  of  the  contour,  c. 

If  A is  taken  to  be  the  planar  area  enclosed  by  the  contour  in  the  xy-plane , for  example, 
one  extracts  an  expression  for  the  spatially  averaged  ''l  z-component  of  vorticlty: 


f.z(t)  ' * U(x,t)  • d?(?) 
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If  a number  of  contours  are  arranged  in  a relatively  fine  mesh  configuration,  the  distri- 
bution of  vvrtlcity  in  the  enclosed  region  can  be  inferred.  1'he  spatially  averaged 
i -component  of  vorticlty  for  the  "i"-th  element  of  the  mesh  can  then  be  determined  from 
the  following  relation: 


C,(t)1  " * ^ <ft  tf(*.t)  • d?(x) 


(41 


1 "1 

Similar  expressions  tor  the  other  vorticlty  components  can  he  obtained  bv  orthogonal 
rotations  of  the  integration  plane. 

The  procedure  involved  in  evaluating  these  integrals  requires  the  accurate  measure- 
ment of  selective  velocity  components  .it  various  spatial  locations.  This  can  t>e  accom- 
plished with  a number  of  velocity  measurement  devices  currently  available.  A sufficient 
number  of  sampling  locations  must  be  employed  along  each  contour  segment  for  an  accurate 
determination  of  the  integral. 

A primary  consideration  for  most  real  flows  of  Interest  is  the  extraction  of  a time- 
varying  mean  value  of  the  fluid  velocity  from  a data  signal  which  also  contains  random 
fluctuations  caused  by  turbulence  and  noise.  In  this  case,  a velocity  component , or 
other  flow  variable  q,  can  be  written  in  the  form 


q(t)  • qftT  ♦ q 1 ( t) 

where  q denotes  the  tlme-varvlng  mean  value  of  the  variable  q,  and  q'  is  a random 
fluctuation  characterized  bv  the  relation 
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The  time-varying  moan  value  of  the  variable  must  then  be  extracted  by  a phase-locked 
ensemble  averaging  technique.  In  the  present  application,  equation  (41  can  be 
rewritten 


VTtT7  ■ - ^ / m?.t)  • d<*(D 


1 
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where  the  overbar  represents  an  ensemble  averaged  time-varying  value  of  the  variable. 


An  evaluation  of  the  entire  tlme-varvlng,  spatial  vorticlty  distribution  involves 
the  determination  of  the  ensemble  average  of  velocity  at  each  measurement  location. 

A sample  of  an  instantaneous  cycle  of  data  is  constrasted  with  its  corresponding 
ensemble  average  tn  Figure  1.  Kach  averaged  data  set  may  require  from  a few  to  several 
hundred  data  "records”  depending  on  the  level  of  turbulence  at  that  location.  The 
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Oscillation  frequency,  f » 8.2  Hi 


Figure  1.  Comparison  on  Instantaneous  Data  Record  with 
Corresponding  Ensemble  Average 
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subsequent  spatial  integration  around  each  contour  requires  that  the  individual  temporal 
ensemble  averages  be  stored  for  later  examination  as  new  measurements  are  taken.  Since 
the  mesh  may  involve  a large  number  of  elements,  each  itself  a contour,  there  is  a pot- 
ential requirement  for  the  handling  of  a massive  amount  of  data,  dictating  a need  for  a 
digital  computer-based  data  acquisition  system. 

To  demonstrate  the  usefulness  of  the  method,  an  experiment  was  constructed  whereby 
an  unsteady,  separated  flow  was  generated  behind  an  oscillating  fence-type  spoiler  on 
one  surface  of  an  airfoil  (NACA  00121  in  incompressible  flow.  The  experiment  was  con- 
ducted in  I'SAF  Academy's  .61  m (2  ft")  x .91  m (3  ftl  subsonic  wind  tunnel.  The  spoiler 
oscillation  frequency  could  be  fixed  at  values  between  0-30  Hz.  Spoiler  mean  height 
and  amplitude  could  also  be  adjusted  for  a given  series  of  measurements . 

A multi-component  constant  temperature,  linearized  hot-wire  anemometry  system  was 
employed  for  velocity  field  measurements  in  the  region  surrounding  the  reverse  flow 
portion  of  the  separation  zone.  An  "x"-geometrv  probe  was  oriented  to  provide  two 
orthogonal  components  of  velocity  in  the  integration  plane.  Analog  output  voltages 
directly  proportional  to  the  instantaneous  velocttv  component  magnitudes  were  then 
sampled  by  a PDF  II/AI  data  acquisition  svstem.  At  each  measurement  location,  data 
records  having  a length  of  one  and  one -ha If  periods  of  oscillation  were  triggered  into 
the  computer  by  an  inductive  pickup  w.  oulstd  o rput  corresponded  to  a prescribed 

value  of  the  spoiler  height  This  process  ha,.  Me  .Meet  of  phase-locking  the  initi- 
ation of  sampling  to  an  event  which  corresponded  o the  periodic  flow  field  drive.  A 
general  schematic  of  the  experimental  apparatus  is  '.'resented  in  Figure  2. 

The  data  reduction  algorithm  involves  conversion  ,f  voltage  to  velocity,  ensemble 
averaging  as  discussed  above,  and  subsequent  storage  of  the  averaged  data  on  disk.  Each 
data  set  was  tagged  with  the  probe  location  and  aerodvnamic  reference  data  (pressure, 
temperature,  free  stream  velocity,  etc. 3,  the  use  of  hot-wire  anemometer  probes  pre- 
sented certain  limitations  on  the  possible  measurement  locations  due  to  the  intermittent 
reversing  nature  of  the  flows.  Data  validation  was,  therefore,  accomplished  using  both 
analog  comparator  circuits  and  a software-verification  scheme  in  addition  to  flow  visual- 
ization photographs.  A complete  discussion  of  these  restrictions  and  the  validation 
scheme  is  presented  in  the  paper. 

Hie  probe  was  moved  and  accurately  positioned  by  a motorized  three-dimensional 
traversing  mechanism  having  a continuous  analog  position  output.  Upon  completion  of  the 
storage  of  each  data  set,  program  control  was  exercised  to  reposition  the  probe  through 
a D/A  converter  voltage  output.  Hits  voltage  was  employed  to  d ive  several  sets  of 
solid-state  relays  in  a power  distribution  circuit.  The  data  acquisition  and  reposit- 
ioning process  was  accomplished  until  the  entire  contour  sequence  was  completed. 
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Finally,  disk  stored  average  data  were  recalled  to  core  memory  for  the  spatial  Inte- 
gration calculations.  Individual  data  points  corresponding  to  common  phase  points  of  the 
oscillation  were  integrated  to  provide  the  circulation  of  various  mesh  elements  at  that 
time.  Repetition  of  the  Integration  procedure  for  other  phase  angles  then  provided  the 
entire  spatial  and  temporal  distributions. 

Results  of  preliminary  measurements  using  this  apparatus  and  technique  are  presented 
for  a number  of  oscillation  frequencies.  A schematic  of  the  rectangular  mesh  config- 
uration Is  presented  In  Figure  3.  Data  from  221  separate  locations  were  sampled  for  the 
31  contour  Integrals.  In  the  examples  provided,  between  21  and  100  records  were  used  for 
the  averages.  Figure  4 displays  the  vortlclty  contained  In  the  contour  labeled  "b"  In 
Figure  3.  Spatial  distribution  for  given  phase  angles  can  then  he  obtained  by  suitable 
cross-plotting  of  these  distributions.  It  is  then  possible  to  track  regions  of  con- 
centrated vortlclty  as  they  move  spatially. 

Tests  currently  being  conducted  employ  a much  finer  mesh  contour  array  and  involve 
a laser  doppler  velocimeter  for  non-intruslve  velocity  measurements  In  the  Interior  of 
the  reverse  flow  separation  vortex.  These  new  measurements  will  also  he  discussed  In 
the  paper. 

Results  obtained  to  date  Indicate  the  capability  of  the  method  to  accurately  determ- 
ine the  vortlclty  distribution  in  a complex,  unsteady  flow.  A major  advantage  of  the 
technique  Is  Its  ability  to  assess  the  average  vortlclty  within  a region  without  affect- 
ing measurements  within  It.  The  spatial  resolution  of  resultant  distributions  is  dict- 
ated by  the  minimum  size  of  the  Individual  contours  which  are  employed.  This  factor  Is, 
In  turn,  often  dictated  ns  much  by  the  capabilities  of  the  available  data  acquisition 
system  ns  by  the  nature  of  the  flow  Itself. 


